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ASEE Report Summer 1982: 

Application of Analytical Redundancy to Failure Detection 

Frank 0. Alexandro, Jr. 

University of Washington 


This project was motivated by the problem of detection of sensor 
failures in a dual tetrad inertial navigation system. The tetrad inertial 
package consists of four laser gyros, three mounted on orthogonal axes and 
the fourth mounted at an angle of approximately 54.7° to the other axes. 

The proposed dual tetrad system consists of two tetrad packages mounted 
some distance apart on an aircraft. Because of the flexibility of the 
structure and the physical separation of the two tetrad units the gyros and 
accelerometers do not measure exactly the same quantities. Ames has 
previously had occasion to study the detection of failures in a single 
tetrad system (Ref. 1) but analytical redundancy techniques were not used 
in this study. The term analytical redundancy is usually taken to mean 
that the dynamic behavior of the system i.e., the differential equation of 
motion are used as part of the failure detection algorithm. 

Several different analytical redundancy systems have been proposed 
(Ref. 2-5) for use in aircraft guidance system. The basic idea is to use 
the information about the system dynamics in such a manner that sensor 
measurements of one quantity can be related via the dynamic equations to 
quantities measured by other sensors. This is in contrast to what is 
referred to as hardware redundancy in which multiple sensors measuring the 
same quantities are used in a "voting" scheme to determine which sensor has 
failed. 

Most analytical redundancy algorithms make use of observers or Kalman 
filters to compare the information measured by different sensors measuring 
different quantities in order to detect failures. During this summer I 
reviewed the available literature on analytical redundancy in an attempt to 
identify those methods which may be useful for failure detection in a dual 
tetrad system. Generally the algorithms described are quite complicated 
and are based on the assumption that the system is linear (or can be 
adequately described by linearized equations). The effects of 
nonlinearity, parameter variation and parameter uncertainty have not been 
extensively considered in published reports. Consequently, it is unclear 
as to whether it is practical to implement an analytical redundancy scheme 
in the dual tetrad system. During the coming academic year and next summer 
I hope to continue to investigate the feasibility of a practical analytical 
redundancy scheme for a dual tetrad system, and I have submitted a proposal 
to NASA for graduate student support on this project. During my stay at 
Ames I have had the benefit of very useful discussions about this project 
with R. Bach, M. Sidar, and G, Xenakis. 
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Does Phosphine Photolysis in Jupiter’s Atmosphere Lead to 
Organic Chemical Evolution? 


Robert Benson 

Nazareth College of Rochester 
Assistant Professor of Chemistry 


My work, this summer, is part of a long term investigation, 

At NASA-Ames, into the possibility of organic chemical evolution 
in Jupiter's atmosphere. These studies have some relevance to 
the origin of life on earth because the earth’s primitive atmos- 
phere may have resembled the present day reducing atmosphere of 
Jupiter. * 

Phosphine, PH 3 , is a constituent of Jupiter's atmosphere. It 
has been detected in the region just above the white ammonia, 
cix’rus clouds that are a dominant feature of our view of Jupiter. 

Light from the sun reaching this region is energetic enough to be 
absorbed by phosphine and cause the photolysis reaction* 

(1) PH 3 + hv--~? PH 2 + H. 

The photons absorbed by phosphine contain more energy than is 
required to break a P-H bond. Most of the excess energy is carried 
away by the H atom as kinetic energy, producing a so called hot 
H atom. It is called this because it has more kinetic energy than 
it would have if it were in thermal equilibrium with the gas. 

Hot atom chemistry, or chemical reactions initiated by hot 
atoms, is a common ocurrance in interstellar regions and the 
upper regions of planetary atmospheres. What I was to investigate 
was the possibility of the hot H atom initiating organic chemical 
reactions in Jupiter's atmosphere. In particular , does the 
abstraction reaction: 

( 2)H + CH4 H 2 + CH3 (feura* producing a methyl 
radical which can go on to produce more complicated organic molecules? 

At the wavelength of light we used 184.9 nm, the hot H atom, 
from phosphine photolysis contains much more energy than is required 
by reaction 2. However the reaction is not efficient and the hot 
H atom can become thermalized by collisions and quickly lose its 
excess energy. Thus it is not at all certain this reaction will 
cfcur^ to any great extent. 

The experiment involved mixing accurately known amounts of 
PHo, CD4 and H 2 and exposing the mixture to the 184.9 nm light from 
a low pressure mercury lamp. CD4 was used in place of CH4 because 
HD will then be formed by reaction 2. Using mass spectrometry to 
measure the amount of HD formed then gives a measure of the 
efficiency of reaction 2 . 

In a simulated Jovian atmosphere consisting of 4 torr PH 3 , 

20 torr CD 4 and 200 torr of H 2 and a time of exposure resulting 
in 15 fo PH3 loss - I measured no excess HD formation over the 
control l A check on the experimental procedure was made by 
exposing a mixture consisting of 4 torr PH3 and 320 torr CD4. 

In such a mixture thermalization of a hot H atom is less likely 
to oour# before reaction with CD 4 . The result was that 0.0088 moles 
of HD were formed per mole of PH 3 decomposed. This finding indicates 
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the experimental procedure is working fine. 

The more realistic Jovian atmosphere results tell us that 
phosphine photolysis and the hot H atom produced is not going to 
lead to the production of methyl radicals by reaction 2. The 
reasoni whyiiis Aprobably because H atoms react very readily with 
PHoi (3) H + PH 3 — > Ho + PHp . This reaction is very close to 
being diffusion controlled and probably siphons off H atoms. 

In summary, my measurement of no abtraction of H from 
methane .under conditions more favorable for this reaction than 
on Jupiter, tells us that phosphine photolysis will not initiate 
organic chemical evolution on Jupiter. 
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WORKLOAD STUDIES IN AFTI F-16 TEST FLIGHTS 


Jan Berkhout, NASA Summer Faculty Fellow 
NASA Dryden Flight Research Center, August, 1982 


During the summer of 1982 I worked closely with several 
groups of NASA and Air Force researchers associated with the. 

AFTI F-16 program. The AFTI F-16 incorporates a number of 
superaugmented and decoupled flight control modes. My aim 
during the summer of 1982 was to familiarize myself with the 
flight test program and data collection procedures so that 1 
could obtain and analyze flight data relevant to "pilot work 
loads" in the superaugmented control modes. 

As used by human factors psychologists, the term pilot 
workload describes an extensive combination of perceptual, 
cognitive and psychomotor elements (Roscoe, et al. , 1978). 

Work overloads are said to occur whenever anticipated performance 
requirements might exceed operator capability. A main assumption 
is that workloads in general ought to be reduced, and that the 
advantages of any augmented control system should consist at 
least partly in reduced workload. 

Attempts have been made to separate motor and cognitive 
workload components (Mornay, 1979). In the operation of high 
performance aircraft a complete separation of the two is not 
possible or desireable. Complex psychomotor functions do 
require physical work, and a reduction in the energy demands 
on the pilot's skeletal musculature often reduces cognitive 
workload as well. 

Psychologists have searched for years for a single parameter 
measure of workload without success. Workload estimates are 
necessarily based on a multidimensional array of information. 
Workload analyses are commonly broken down into four broad 
categories as follows: 

a. Primary task analysis (Chiles and Alluisi, 1979) 

b. ' Secondary (probe) analysis (Ogden et al, 1979) 

c. Physiological monitoring (Wierwille, 1979) 

d. Subjective reports (Williges and Wierwille, 1979) 

These four categories often generate independent and 
contradictory estimates which are not easily reconciled 
(Hart et al, 1982). Workload seems to be a genuinely multi- 
dimensional consruct, and no one category of measures can be 
adequately interpreted without information from the others 
as well. 
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Work load measurements In all four domains are being Incorporated 
into the current AFTI F-16 flight test program . (Crombie and Frazier, 
1982). The main interest is, of course, comparing work loads imposed 
on the pilot when the unconventional control modes are used, and the 
transition phases from one mode to another. Further on in the test 
program a number of automated attack and maneuvering systems will be 
made operational and studied, Including voice displays and controls, 
and the entire rational behind these systems is in fact a reduction 
of pilot workload. 

My main associates during the summer of 1982 were the AFFTC 
personnel coordinating Human Factors evaluations in the test flight 
■program, particularly Maj. S, Gray and Capt. R. Crombie. P. Sorokowski 
was my chief contact in the Joint Task Force (JTF) operations office, 
making flight data available to me, and briefing me on mission 
scenarios and data formats. The NASA-Dryden researchers in the 
Dynamics and Control Branch were valuable sources of information 
on flight dynamics and control theory, particularly P. Carr, and 
from AFFTC, B. Schofield and T. Twisdale, Finally, on several 
occasions I was able to meet with Dr. S. Hart at NASA-Ames for 
discussions on more general applications of work load measurements 
in aircraft operations. 

Primary task analysis . Primary operator outputs in the AFTI 
F-16 consist of a steady stream of motor activity into four primary 
flight control axes (Sr; Sp; P and T) , the T control being a throttle 
twist grip not present on conventional aircraft. In some cases 
the operation of the primary flight controls can be described in terras 
of continuous operator transfer functions (Heffley, 1982), for which 
a finite number of gain and lifhase lag values can be specified. 

However, in many cases pilot output cannot be adequately specified in 
terms of continuous transfer functions, and is better handled in terms 
of Integrated Sequences of Activity (ISAs, Berkhout, 1981), or even 
as conditioned responses. In such cases the appropriate parameters 
to monitor would be response time histograms (Berkhout, 1982), or 
error values relative to information content ( Shingled ecker et al, 

1980) . 

When closed-loop tasks, such as tracking a known target, are 
performed an RMS error value can usually be calculated for the 
concurrent operator transfer functions (Poulton, 1974). For open 
loop tasks, where target position or the value of dlsequilibrating 
inputs are not known, error functions cannot be quantified. However, 
it is possible to evaluate performance and workload in such cases 
using spectral density analysis of output in each control channel 
(Sheridan, 1974), evidence of single channel time sharing (Kinsbourne, 

1981) , and features of control outputs related to queueing theory 
(White, 1975) . Applications of these techniques to high performance 
aircraft studies are discussed in Zacharias and Levinson (1979) and 
Schofield et al, 1982). 
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My analysis of AFTI F-16 flight data for the Initial flights 
in the program Is based on the manual reduction of control force 
nnd position values derived from strip charts of telemetered 
parameters recorded in real time during each flight (Table 1). 

Epoch boundaries are determined from other flight parameters 
on these strip charts as shown in Table 1, The value of 
subsequent work load estimates depends on the correct determination 
of epoch boundaries, which should indicate such things as the 
initiation of specific maneuvers, changes in control modes, and 
arbitrary events associated with maneuvers such as reaching or 
exceeding certain angles of attack. ■ 

Once epoch boundaries have been determined, the strip chart 
■data is examined to establish values of control activity in each 
discrete axis. This can be done by noting any above-threshold 
control events, and tabulating reversals of the primary controls. 

Once this has been done, the activity in all control axes can 
be collated and examined for evidence of parallel and sequential 
control outputs. From such data, evidence can assembled concerning 
the operators control strategy expressed in terms of cycle time, 
sequence patterns, switching rates and dwell times, and from these 
-data inferences can be drwn concerning the relative work effort 
across epochs. 

In the three flights flown up to August 24, 1982, only the 
standard normal control mode has been used. Examples of strip 
chart data are given in Figures 1-5, and these indicate that 
full time sharing of the three primary control axes occured, 
as would be expected with normal control configuration. 

Secondary task analysis . While my main effort this summer has 
been to set up analysis procedures for the primary control data. 
Information from the other three categories of work load measures 
will eventually be available as the flight test program proceeds, 
and will be collated with the primary control output data. 

Traditional secondary task analysis includes such "probe" 
tasks as voice shadowing, reaction time to extraneous stimuli, and 
mental arithmetic. These are presumed to reflect reserve processing 
capacity when performed in conjunction with primary control tasks. 
(Note that only the instructions to the subjects provides the dist- 
inction between primary and secondary!)" To be effective, the probe 
tasks ought to impose a constant load at a fixed pace, and meet 
several other constraints as well (Brown, 1978), For these reasons 
probe tasks as such have not been included in AFTI F-16 test flights. 
However, it is possible to regard some forms of voice communications 
and flight controller instructions as secondary tasks with a precisely 
defineable information content (Shingledecker et al, 1980). When 
such communications occur spontaneously during complex test flight 
scenarios response time and accuracy can be scored, and the impact 
on primary task functions can be quantified. 


i 


ORIGINAL. PAGE !3 

“ — /MtAl ITV 


Physiological monitoring . Most physiological subsystems 
show reliable changes with increased cognitive and motor work loads. 
Under high work load and stress conditions cardiac functions, 
energy metabolism, blood chemistry and central nervous activity, 
are all rather highly correlated (Berkltout, 1970). Cardiac 
functions will be monitored on some future AFTI F-16 test flights 
(Reader, 1981) . Basic rate, rate variability and qrs / t wave 
phase relationships will be recorded. These physiological 
parameters will be available for correlation with the estimates 
of work load from the other three categories of data, and will 
provide a useful interpretive aid. 

Subjective workload estimates . There is a large literature 
of special purpose workload estimation instruments designed for 
utie with high performance aircraft operations, which combine 
multi-attribute theory and scaling techniques to produce equal- 
interval response matrices for grouped pilot responses (Donnell 
•and O'Conner, 1978; Reid et aT., 1981^ Williams, 1982). These 
response measuring instruments are in theory to the traditional 
ordinal scales of aircraft handling qualities such as the Cooper- 
Harper (1969), a t»en rank rating list. The best of these multi- 
attribute instruments unfortunately require a large amount 
of interviewer-pilot contact time, 'and a pool of group data. Both 
these requirements are unrealistic in terms of the AFTI flight test 
program constraints. The emphasis in AFTI debriefing will be on 
segregating pilot workload estimates for single mission phases 
(Greene et al, 1981) rather than obtaining load estimates for 
• larger mission segments. The methodological emphasis will be on 
obtaining a true, discrete boundary for the desired event or 
control condition, which may be quite brief. 

Continuation studies . As more test flights are completed, 

Major Gray and Mr. Sorokowski have agreed to send copies of the 
appropriate strip charts containing pilot control outputs to me 
at the Human Factors Laboratory at the University of South Dakota. 
Further analysis of the strip chart data will continue during 
the academic year 1982-1983 along the lines indicated above. While 
I will be concentrating on primary task data, information concerning 
work load estimates for the same flight epochs but derived from 
the other three categories will be made available for correlation 
with the primary output data. The primary data analysis is 
generally based on an open -loop interpretation of the test flight 
scenario. Closed-loop data may become available later in the 
flight program when tracking scores for 1IQDT or "aggressive tracking" 
become available, and workload scores can be collated with tracking 
values. In particular, the early flight data must really be 
considered baseline data, since the augmented control modes and the 
fourth control axis have not yet been used. Operator output bandwidth, 
channel capacity and queueing of output elements will be of much 
more interest in the four axis augmented modes. 
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♦ 

’ Apart from the further processing of strip chart data 
several other procedures for analyzing AFTI F-16 workload 
levels are possible. Power spectral density and cross-spectral 
transformations of the motor output time series would permit 
a more precise look at bandwidth and shared power aspects of 
multi-axis tracking. Once available, these transformations 
can also be used for calculating event-linked postgrade and 
retrograde averaging of control operations. In such a data 
reduction, the "events” would be critical points in the 
flight scenario, and the averaged responses' would show such 
things as the time needed to damp out responses in particular 
frequency bands. 

Once the magnitude and quality of data generated by the 
AFTI F-16 test flight program is known in more detail, I will 
present a proposal for NASA support of further work, along , 
the lines suggested in the paragraph above, to be done at 
the Human Factors Laboratory at the University of South 
Dakota, and involving approximately 12 man-months of data 
..analysis effort by senior graduate students in the doctoral 
program in human factors, » 
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Infrared Observations of M17 


Elizabeth Brackmann 

State University of Hew York at Stony Brook 


There is increasing observational evidence that massive stars 
tend to form near the edges of dense clouds of gas 
and dust in the galaxy, leading to the formation of "blisters" of 
ionized hydrogen at the surface of the clouds. This summer I 
have studied infrared emission from one of the most energetic of 
those star forming complexes, M17. Attention has been focused on 
the properties of the emission at 20, 50, and 100 microns, on the 
temperature distribution of the IR emitting dust, and on the role 
of several stellar-like point sources recently discovered in t^*e 
near-IR by Ames observers in the overall energy balance of the 
region. The results of this analysis have been combined with the 
work of past investigators who have studied other aspects of 
M 17 — its stellar populations, radio emission from diffuse and 
compact regions of ionized hydrogen, millimeter wavelength 
observations of molecules deep in the cloud core and IR emission 
at other wavelengths and resolutions-- to refine 5 our view of 
ongoing star formation in this region. A paper is now being 
prepared on this work in collaboration with Mike Werner of the 
Space Sciences Division. 

The initial analysis was done using the high resolution (30 
arc seconds) data at 50 and 100 microns taken with the Kuiper 
Airborne Observatory. Dust temperatures were derived by fitting 
a black body temperature curve to the 50 and 100 micron emission 
at each measured point and then correcting for a dust absorbing 
efficiency which goes as 1 /k. The temperature decreases steeply 
from east to west. Comparison of the 100 micron brightness with 
that estimated for a blackbody at the observed 50-100 micron 
color temperature indicates that the opacity increases steeply to 
the west. In addition, the 50 and 100 micron emission maxima are 
displaced from one another in the direction of the steepening 
density gradient. The most likely interpretation of these results 
is that a high extinction cloud is heated from its eastern side 
by young stars located near its edge. The longest wavelengths 
penetrate the furthest into the cloud with the least ability to 
heat dust, hence the warmer (50 micron) dust lies closer to the 
heating source than the cooler (100 micron) dust. 

This model is an extension of the picture proposed by other 
observers. The 6 cm radio emission from ionized hydrogen 
adjacent to the far-infrared peaks has been interpreted as the 
front between the HI I region and the molecular cloud seen edge 
on. The present results depart from previous work in showing 
that the peak IR luminosity in the molecular cloud is too bright 
and compact to be heated by the star clusters to the east which 
• are thought to power the H 1 1 region. My hypothesis is that there 
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is a deeply embedded unobservable star cluster closer to the IR 
peak providing local heating. The newly discovered near-IR point 
sources lie in the predicted location for the new heating source. 
Examination of the spectrophotometry of the 14 near-IR sources 
shows that they are very luminous but cool { 100°K) objects, most 
likely massive young stars ( 30,Q00°K) throwing off dust shells 
heated to 100* K. The heated shells may be the only evidence of 
other nearby young stars without heated shells, but highly 
obscured by foreground dust, whose unbounded radiation is 
powering the local area. Work is now in progress on the other 
possible interpretations of the near-IR point sources. 

One of my main motivations for coming to Ames this summer 
was to become more familiar with the techniques of IR 
instrumentation and data handling. I have had several useful 1 
demonstration-discussions with my NASA colleague Tom Roellig on 
the state-of-the-art research being done in his lab (NASA-Ames 
Space Sciences Division) on low temperature high sensitivity IR 
detectors. During the summer I also spent a substantial fraction 
of my time reducing the data for a very high resolution (5 arc 
second) 20 micron map of Ml 7 . The data was taken using a 
beam-switching mode in which a position is sampled at the same 
time as adjacent positions a few arc seconds on cither side. 

When the adjacent measurements are of blank sky they provide an 
easy means of removing atmospheric contributions from the signal 
at the central point. In this case adjacent measurements were 
taken within the Ir emission region and a scheme had to be 
devised to deconvolve true signal from the relative measurements 
at each point. The final 20 micron map shows a number of peaks in 
addition to an extended background which is co-extensi vc with 6 
cm radiation from the ionized hydrogen at the edge of the cloud. 
The 20 micron peaks coincide with compact HI I regions and with 
some of the near__IR point sources and have been useful in 
determining the nature of the new embedded cluster. 
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1 . ANALYSIS OF TWO-WAY CONTINGENCY TABLE BY ROCK AND VEGETATION TYPES 

2. MODIFIED MAXIMUM LIKELIHOOD CLASSIFICATION PROCEDURES AND THEIR 
APPLICATION TO PATTERN RECOGNITION IN REMOTE SENSING 


Hubert J. Chen 

Department of Statistics and Computer Science 
University of Georgia 

1982 ASEE Summer Research Fellow 
Stanford University and NASA/ Ames Research Center 


Dr. David. Mouat at the Technology Applications Branch is currently conducting 
a- research project in geology on determining the relation between rock 
and vegetation types. The sampling survey and statistical analysis are to 
be used as basic tools for the study. It is suggested that a minimum 
sample of 250 items for the ten by five contingency table be collected by 
a random process. It is also suggested that a large sample be considered 
if more budget funds and time are allowed. 

Due to possible empty cells in the contingency table, the actual statistical 
analysis could be more complicated. Empty cells occur in two situations: (1) 
when the sample size is not large enough (small sample effect), and (2) when 
the cell combinations are impossible (with probability zero cells) . A reso- 
lution for the former case is to add a small number, say 0.5 or I to every 
cell and then the usual chi-square test for independence (or equivalently, no 
interaction) could be employed [J.) . In addition, a log linear model via 
logorithraic transformation of the expected frequency may be used to convert 
this table into a two-way ANOV A type problem. Thus, Tukey's one-degree-of- 
freedom test for interaction is possible and furthermore, test and compari- 
sons among row/column effects can be made. For the latter case the with- 
probabillty-zero cells cannot be adjusted or estimated; a special method 
through iterations must be used to derive the expected frequencies for the 
nonempty cells. Then, a chi-square test for the quasi- independence can be 
performed [1] . Further analysis can also be made by a log transformation for 
all nonempty cells. In this situation an ANOVA test is possible using re- 
gression model via the use of SAS statistical package [5] in which sum of 
squares type I, II, III and IV should be read and interpreted with caution. 

The second problem I have been working on is pattern recognition in remote 
sensing via a satellite, which was hinted by Dr. Don Card of the Technology 
Applications Branch, NASA/ Ames Research Center. This subject interests me a 
lot. A modified maximum likelihood classifier which should have the potential 
of increasing the accuracy of classification is studied. 

At the present time statistical model used to classify a picture element 
(pixel) on earth into an appropriate category of land cover and land use is 
a multivariate normal probability distribution on the brightness values or 
color-coded values (X^, X 2 , •••, X^) collected and digitized by multispectral 

scanner system (MSS) installed in a land satellite [3, 6], The domain of 
the model which is assumed to be in the entire multidimensional space may 
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not reflect the practical situations. Instead, it is more likely to be 
restricted in certain finite intervals say (a,b). The brightness 
values from a spectral band, for example, are confined in the range of 
(0, 127) in an MSS of Landsat. Based on this reason a doubly (two- 
sided) truncated multivariate normal probability model is assumed on the 
vector (X^, Xg, . .., X^) ; and the conjecture is that based on the trun- 
cated model the accuracy of classification of a pixel into a class could 
t>& increased using maximum likelihood classification principal. 

Since the true parameters of the probability model are generally unknown, 
one has to estimate them using training samples of known features on the 
ground. To do this the maximum likelihood method or moment method can 
be used [2, 4'J. As a result, a system of simultaneous p(p+3)/2 nonlinear 
equations must be solved by high speed computers. A generalization of 
the model to several classes with different mean vectors but with the 
same covariance structure will be investigated in the future. Comparisons 
between nontruncated and truncated probability models will be made 
accordingly in terms of the chance of misclassif ication. 
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Facilities for Simulation 
of 

OTV Aerothermodynamics 

Dick Desautel 
Associate Professor, 

Mechanical Engineering 
San Jose State University 

My summer project has consisted of developing an understanding of the 
problem posed of experimentally simulating orbital transfer vehicle (OTV) 
aerothermodynamics. The proposed OTV mission involves a unique flight regime 
for the dominant aerodynamics and heat transfer. This flight regime is one of 
hypersonic (18 to 37 for M) , low density (Kn =* 0.001 to 0.3), and intermediate 
Reynolds Nunber (10 to 100,000) flow. The flow is dominated by real gas effects and 
and nonequilibrium conditions due to the high total enthalpy. This regiinf, has 
been encountered on previous missions , but has never been the sole regime in 
which the vehicle experiences the significant atmospheric effects. The summer T s 
project addressed the question: What is the experimental aerothermodynandc 
.simulation problem posed by this mission and how should we go about solMngit? 

It will be. recognized that the state of knowledge concerning gas kinetics, gas 
dynamics, and gas chemistry is far from adequate in the hypersonic, low density 
regine, and hence there is a strong motivation for achieving experimental, 
simulation of the aerothermodynamics. 

During work on the project, I reviewed the literature, both on fundamentals 
and on applications, of impulse-type flow facilities with emphasis on shock tube 
derived facilities. I also interacted strongly with Dr. Chule Park to gain 
insight into particular advantages /di advantages of each device, and pragmatic 
performance capabilities. At the same time, I studied the basic concepts of 
similitude as applied to the hypersonic and low density (rarefied! gas regimes. 
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A further background investigation considered the rationale for OTV flight regime 
and its consequent free-stream flow parameters. 

All this background and developing familiarity with the problem was synthesized 
and a simulation strategy evolved. The project's conclusions were primarily that 
separate simulations of flow field (aerodynamics) and the relevant gas kinetics/ 
chemistry were needed, the -former involving first simple glow-discharge diagnostics 
in a simple shock tunnel and the latter involving long-term commitment to proper 
modification of existing shock tube facilities to produce the required gas purity 
and total temperature. 
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"Model Reduction I loti iatJs for Manuaver Autopilot Desirin'* 
Kennel h bun i pace 

IiuJi ana Uni ver si tv Purdue University at Indianapolis 
IT o f ussor o f Ei i ci :i r i ©or i. nq 


Tl x.: Pi ol d emu 

Linear dynamic models of aircraft and other complex systems 
are often represented by a hi qlv-or der differential equation. Such 
models are ususally derived from complicated analytical design 
p r or: eel ur oh and involve predicted values of numerous system 
parameters. Tin* complexity of such models provide a serious 
handicap to under standi nq system behavior. Similarly, the 
ouinu I at j vo of Poet 01 uncertainties in predicting parameter values 
lead to gross inaccuracy and misleading characteristics in the 
high order models. In addition., 3 ow-orcler models are reported to 


be i mpor tun l In the design of r ubust controls* 

llio objectives of this project is to develop useable and 
reliable procedures for extracting from a high “-order model those 
1 owei t.n dei models wIM. c.li r et ain the ©ssetial dynamics of the real 
system and provide desi r oabl © insight, or design properties. Such 
procedures would be of current, value in the design of control laws 
foi the ilanuevri Autopilot being used in. the Hi •"•Mat aircraft, 
they will be of similar value in i mp3 ©mentation of control laws 
for future application of the Man never Autopilot concept. 


St. a Lush 

During tbi 


summer, .1 have initiated a literature review., 


applied each of four different methods to a simple example . system 
to illustrate the f unot.J uni rig of the method. In addition, I have 
spent Lime learning to use the Cyber Inter active Facility and the 
ORmCI.S subroutines in applying the methods studied and in 
analysing methods still to be considered. 


Li ter a Lur e Rev i ew ii 

At present.., I have obtained ten technical reports dealing 
with different methods of model reduction. These reports include 
references to additional reports. Thus., there seems to have been 
a significant amount of study on this subject. I expect to 
continue m/ review of existing methods. 

Throe of the methods presented < FERNANDO , MOORE., and PER ME BO ) 
are ape-lied to systems which have previously been "balanced". The 
concept of "balancing" is presented by Moore. His report was one 
of tiro first obtained and I have spent much of the summer 
developing the a! a or i thm and associated computer program to 
pei form tire balancing procedure. The "balancing" procedure may be 
iniri hpical 1 y valuable as it is claimed to lead to a robust form 
of control 1 aw. 


Iloor©* « report contains 31 references, some of which should 
bo rovi owed for additional information. The report also contains 
riuiiiwi ou5 proofs to support the development of the concept and 
sovc-i al examples of its: application. the Moore method of 
r eduction is one of the four which I have illustrated by example- 
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Imt! to ( In'* 1 tw lit:* amount of time Bpnrvfc e*«tr at t inq the 
"bui , tr it* i mi" uluoi iUun and dove! oping the computer program to 
:uup3 omen t; l ho al»iurlt I itn 1 have not yut reviewed tiio oilier two 
mol 1 uulo applied to lhu "balanced" system (FERNANDO, PERNKI30) . 

o recent i- opart (GUTMAN) describes c\ rather simple method 
In i sod on di Hei enl a al i on of pol vnomial s„ it was selected for 
t oview lnn.(.uiwi;> j i was* si mpl e .and I needed material For my mid*' 
sumitiei pi w&ent.ut i on. It is one of the four which I have 
i J 1 usi rested hv ». xamp.1 o. In i l o present Form it applies only tu 
si null" i nput. simile -output, systems, An extension to the mul ti ~ 
variable (asi- should bo t onsJ rler ud. 

In s.« i i.vi;iil ctnd- (WAI.KLR) fui lif-RF, a met hud described as 
"i era due i * ■ l 1 <■■?! i i ,i tin" is used ns port of the process of ol.rl.aini nci a 
i tnluciMl model t.d the combi nod svsiem composed pf the aim aft and 
the pi J tit. . in t.l'rj r> 'd iulv, the alnl urdei model is reduced to a 

'Order model .. It ion’’ t clear whether this method was 
or ivii noted by Walker of cl eve 1 oped form another- report (SKELTON) . 
nuilo it lot of 1.1 mo was spent stuchcimi the Wallen report and the 
"land due retention" unit. hod, 'this met had is one of tire four which 
] have illustrated by tsx anip 1 c*» • .1.1 liar- the undecir cable? property 

of I nt. i ovlur :i nt i spur i pus hi qh f r equerry zeroes into the recluced 
c.-\*totoii» l.rasf er Function, ,l. have* not reviewed the Skelton rep or t >M 
but expeul that, such a review w:i j 1 clarify the origin of the 
"residue r etent i. on " method,, 

The Four ill method Hlubti ated by example was called "State 
Variable Observation" and is si mol y a stat e~-vai- i ab 1 e 
representation of the "Domi nar.t-Root Concept" presented in 
eleentary control textbooks <e. tj. DORF) „ 

^ _ Ihe other reports (FI -ATTAR, KOKOTOV It:, KRXSHNAMtJRTHY, and 
SNAKED) 9 which have not been reviewed;, are listed in the 
hi I f.’l :i ecu 'upi > v. 

ftio La amp 1 e Cyst urns 

Tlie example system „ used to illustrate eatt.li of the methods, 

i s s 

I i n > *-0 . 4B (b KG ) ( s i-ki ) / < i- 1 ) c s 1*2 ) < &+6 ) 
fliis system lias an i. mpul. s,e response* of s / y. 

■#sC 6 ) ~ ~ &*** *+■ 0 ,ai<f£~ s ^ 

and a slop rospoi.so of s , 

^CsPl^) -/" - /# 5 f¥ £ 4 ~ <3 * && £ ~Xjc w o,oo H( 2 T 

I he Cvstorn was chosen because, clear Iv, the response clue to the 
eigenvalue at 6 contribute?: little to the* response and the system 
diu-s Ivor a quod lower • order model,, Hence it should bo effective 
in displaying the pur for manor of various methods. A set of view-* 
i.ii wphs was pr f.-pat ed womiiai inn state var '.table models,, 
ti a i is f oi t uiiu t. i oi i model. s„ i iiifiuJ. so i '"c-'Sjiousos, and step reBponsres of 
this system usirn) the classical "Dominant Root." method (DORP) and 
each of the fcuu methods described above. 

Duo lo Fjpaco 1 i mi tat ions, only the stop responses are 
pi t .-rt -II I Od in sumuiai ys 

A„ i-'orni nan t - R( Kit : »• . 

a J-aa. M - / - /, s'Ve “* *4. o e ^ 

I*„ .it air? Vari abl r ni.u”£M vation;i 

$m.o, o= /- /.s-ve - ^ + o,s c te~ ) -* 
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C„ Residue Rolenl :i tins , -jt 

#«.<*> = 1-J<r3e. 4-o,s-j& 

D« Differ Dirt i ati on of -Polynomial©* 


F“„ Moan"? methods 


i-» », m. vir p i v» » I v* •*» » i hjiii* i.' «• »• u j 

4 *>u)~ }-hyte~ ,,X!rt -t-t.zie 


- 3 , 1 ')* 


4^U)~0,73 -M l /£- , ' ox± +0.3i>Zd 


rlfiox 


Si at o Vur i ab 1 e Recover y s 

Most of the methods of mud el reduction involve similarity 
transformation of the system model . Tims,, tin- reduced model is 
o>; pt <• nsiietJ :tn terms of stake vuriabl tvn which differ -from those of 
ihe ori qinul system model . If, as is> usual 3 y the cane, the 
vai ■ i c«j 1 es> in the cMiqinal model at e physical Vctriulili'S of real 
.in tea* not,., thou a method must be derived to recover the original 
p tutus from the roilucod state**. The si mi larity ft an •“•formation 
cannot Ik-; reversed because of the reduction. hs an initial 
at tempt at sol vino thi e» problem, I developed u proqram which 
reduce© sy stains winch have; been "balanced " an if then recovers a 
1 esat ■ • squared esti mute of the ciriqinal states. While the method 
wort.:* malhemal: i cal 1 y , it does nut sworn to produce physically 
reasonable rusuJ tu« Further study is needed. 

Ui hi / otiraphys 


DfiRF i R:i char d 
J 900 


hoi f .1 "Modern Control Systems" , Addi son -Wesl &y 9 


EH -ATTARs R. A. FI Attar and H. Vi dyasaqar , "Order Reduction by 

3 1 and 3 Norm Minimisation",, 1EEI1: Transactions on Automat i c 
Control,, hoj. . AC--2Sj, Wo. 4,, Auq, 78 , pp 73 1-/33 


FERIli iML'Q; 


K„ Vi neons; a Fernanda and IT,, Nicholson 


Pertui hati ona.'t Mode,! Reduct i 
Transactions on Automatic Ci: 
72 >i PP 466-468 


'Si naul ar 


on of Balanced Systems",, JEFF 
:«ntrol , Val . AC-27, Na. 2, Apr. 


SUTI li-iNs Per -81 of Gutman ,, Curl Frederick Mariner f el t , and Per 

Mol an rlei • ,, "Contributions to the Model Reduction Problem"* 
TERR Transact i ons on Automatic Control » Val . AC-27 s Mo. 2, 
Apr . 82,, pp 404-455 


KOKH TOv'i i'.t P» V., Kakatovic, R. It. CF Halley, and P„ Sanrmt i , 

"Si nctular Per tuc Uat,i oris and Model Reduction in Control Theory 
An Ov«i *• i no" . Automatic©, Vai . 12, 1976, pp 123—132 


KR 1 8HMAI1URTHY s V. Kri rhnamurthy and V. Sasha;’ i , "Model Reduction 
us:i nq the Ruu I: I r Stability Criterion",, IEEE Tr an sac t i on s on 
Automft fc i c Control ? Vo 3 . AC-23,, No,, 4., Auq, 78, pp 729-731 

IIGORFs Bruce c„ Mour e,, "Principal Component Analysis in Linear 
Systems* Control lability, Observability, and Model 
■Re duct i on 1 ,, IEEE Transactions on Automatic Control , Val . 

AC 26, Mo. i „ Feb. 81., p|.i382-3B7 


Pt:lilli:.nci! Lars Ivrnclio and luunard II. aily«r«an, 

Via Balanced Stakes Bpacn tenrKpkntwna . ifefib l> any.'U.»un- 


tin Automat, i i' Control, Vol » AU' v:? » 
302 ** 3(37 


No. 


Apr * (32. , pp 


S! lAKKlis U, Shaken! and M. Kart.miii.ts. "The Use of ^ros ant 

Zern~M recti on© in Modal -RetJur.t. i wn " , International t uiu naJ < f 
Control Vnl « 23, No. l 9 PP US-UI1 

Ol-Tt TONS Robert. K* Skel ton « "'Cost. Denomposi t i on o\ Linen 6 /stuii » 
j ™U, Appl j cation** l-n lludni Roductlor.% _Intorna«onal Journal 
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Microbial and Nutrient Analyses in Controlled Environmental 
Support Systems for Plants 

Diane Dudzinski 
Professor of Biology 
College of Santa Fe 

if higher plants are to be used to supply food, oxygen and potable 
water as life support for man during large scale manned space missions, then 
it seems reasonable to assume that bacterial accumulations will occur within 
such systems. These potential contaminants can enter via the air, the nutrient 
medium, man and plants. They may have harmful and/or beneficial effects depending 
upon the species. Therefore, it may bo necessary to control microbial popula- 
tions within such systems. In addition If nutrient solutions are to be recycled 
within such systems, it may be necessary to know at what rate specific nutrients 
are depleted from the system for different plants. 

Steve Schwartzkopf of NASA/AMES has developed a computer monitored closed 
plant chamber which can precisely monitor and control 3 number of physical and 
chemical gas parameters. I came to NASA/AMES this summer to work with Steve 
Schwartzkopf and Robert McElroy In the CELSS program, and to investigate the 
above mentioned problems. 

Experiments run this summer include the growth of 'Grand Rapids' lettuce 
from seedling to harvest in an open environmental plant chamber over two, four 
week periods. During these experiments bacterial and nutrient samples were periodi- 
cally collected for analyses. During the first four week period a UV lamp was 
continuously "on", and during the second four week period a UV lamp was continuous- 
ly "off". In another experiment involving the closed plant chamber, a tomato 
plant was placed Into the chamber and samples were periodically taken for bacterial 
and nutrient analyses. 

All microbiological media, supplies, incubators, Quebec colony counters, 
plant growth chambers as well as atomic absorption spectrophotometer and reagents 
for nutrient analyses were supplied by NASA. 

Bacterial findings to date show that UV light emitting wavelengths of 
245 nm can reduce bacterial populations from 60,000 cells per ml. to less than 200 
cells per ml. within 24 hrs.per 100 liters of nutrient solution having a flow 
rate of 20 l . per min. (see FIGURE 1). Bacterial species identified within the 
chambers include pathogenic Staphylococci spp., Erwinia sp. ; opportunistic microbes 
such as 1 Pseudomonas spp., Chromobacter i urn sp . , and other microbes such as Lactoba- 
cillus sp., Acinetobacter spp., Sphearoti ius natans , Leptothrix sp., Serratia 
liquefaciens and Aical igenes sp. . Several of the species identified are involved 
in the production of ethylene and other gases as well as toxins and could be con- 
sidered phytopathogens. Nutrient analyses to date are incomplete, but do suggest 
that metals such as iron are rapidly depleted from the nutrient medium and may 
oe complexed as precipitates onto the sheaths of Sphaerotilus natans and/or Lepto- 
thrix sp. Completion of the analyses from these experiments will be undertaken at 
the College of Santa Fe during the year. 

If funding is approved for a second summer as a NASA/ASEE, I would like to 
examine the same parameters under intermittent UV light exposure which may be 
more energy efficient and just as effective in controlling microbial populations- 
In addition I would like to examine interactions between higher plants, algae and 
microorganisms within a closed environmental chamber. 




DAYS OF EXPERIMENT (7/6/82-7/23/82) 
CLOSED AEROPON1C PLANT CHAMBER WITHOUT UV LIGHT 


FIGURE 1. Microbial cell numbers In controlled environment plant chambers 
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THE PAY-OFF OF HIGH LIFT 

Donald L. Elder 

San Jose State University 
Lecturer in Aeronautics 


This project wa.s to examine the potential benefits of 
further developments in the area of conventional mechanical 
high lift aerodynamics for subsonic aircraft. Power augmented 
* lift was excluded. This topic was suggested by NASA colleagues 
V. Cor si glia and Larry Olson. 

The approach taken was to review the literature for the 
state of the art in aerodynamic high lift technology, for 
aircraft design optimization with respect to high lift and Cor 
the economics of aircraft operations with respect to the 
significance of high lift capabilities and then to draw some 
inferences as to the improvabill ty potential and the cost/benefit 
aspects of advanced high lift performance. 

NASA/AMES is a good location to carry out this investigation 
because of the key research being done here along with the excellent 
library/search facilities. Theory, powerful computerized analytical 
methods and outstanding research test facilities were found to be 
producing exciting advances in high lift. NASA sponsored studies 
by aircraft industry contractors nroduce excellent independent 
assessments of the state of the high lift art and its 

appropriate applications to the aviation field. Included in this 
review were studies by NASA and contractors of High Performance 
Single Engine Airplanes, Applications of Advanced Technologies to 
Small Short Haul Aircraft, Applications of High Technologies to 
Commercial Advanced Designs and High Lift Concepts Selected for 
the Energy Efficient Transport Program. 


Aerodynamic theory suggests that today’s maximum lift and 
lift— to — drag ratios could ideally be doubled or more. Nearly 
half of what can be demonstrated on an experimental basis for 
two-dimensional flow is lost when applied to a complete aircraft. 
Nevertheless, great advances are being made in increasing two- 
dimensional high lift. The heretofore unusable region beyond 
the onset of wing flow separation may now become part of the 
designer’s province for tailoring an optimal aircraft design. 
Normally unacceptable flying qualities near, at and beyond the 
stall have resulted in conservative Federal regulations which 
effectively prohibit utilization of from thirty to forty per cent 
of the available well behaved upper portion of the presently 
demonstrated high lift capability. 


A spectrum of aircraft types and operations were examined 
which included single engine high performance personal/business 
aircraft, small multi-engine executive/busincss aircraft, medium 
haul ^0-l£0 passenger aircraft, short haul 3G--60 
and long haul major air carriers. The . high lift 
tune were found to be different, especially 


passenger 
needs for 
in emphasis. 


aircraft 

each 
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The economic factors are initial cost, direct operating cost, 
total cost per seat mile, productivity in seat miles per hour per 
dollar, utilization potential in revenue miles per day, break even 
load factor and the particular market being addressed. Implicit 
herein is the extremely powerful factor of ever Increasing prices 
of petroleum fuels for which in the foreseeable future there is 
no alternative. 

The typical flight profile may be divided Into several essent- 
ially Independent segments with respect to the aircraft’ 3 high 
lift configurations: take-off ground roll, take-off climb to clear 
t an obstacle, second segment take-off climb with one engine Inoper- 
ative, clean climb to cruising altitude, cruising flight, descent, 
landing approach and touchdown, landing ground roll and alternat- 
ively, climb following, an aborted landing. Each of the foregoing, 
under some circumstances, can be a limiting condition with respect 
to aircraft utility depending on factors such as field length, alt- 
itude, temperature, aircraft design characteristics such as power 
to weight ratio and. the wing design, and, not the least, the 
constraining influence of Federal Aviation Regulations for insuring 
safe flight. 

Aircraft lift and drag are proportional to the wing size and 
the airspeed squared. The constants of proportionality for lift 
and drag are and C^, varying with the angular attitude of the 
flight. Maximum lift at Crmax would provide the lowest landing 
or take-off speeds and distances. Climbing performance is most 
strongly related to the lift to drag ratio L/D. Best range is 
obtained when the product of the speed times the lift to drag ratio 
is a maximum (VxL/d)max (for jets). Sophisticated aircraft designs 
may provide a different high lift configuration for each segment 
of flight. 

Pay-off potentials for a five per cent improvement injhigh 
lift for a flight segment limited jet transport can typically be 
a twenty per cent increase in payload at take-off, a forty-five 
por cent increase in useful load in second segment climb and a . 
fifty per cent increase in payload when landing weight limited. 
Spectacularl however, it was found that whereas most aircraft 
specifications/designs call for maximal high lift capabilities, 
most aircraft operations are not constrained by high lift 
limitations. It would appear that If all existing aircraft ’were 
to suffer a five per cent reduction in their high lift capabilities 
tiie economic impact might go almost unnoticed. What then is driving 
high lift research and development? Among other factors aircraft 
are being inefficiently used on short hauls for which their con- 
figurations are not appropriate and the airlines are losing money, 
but deregulation will result In demands for the service being 
met by technologi cally and economically appropriate new aircraft. 

Low initial cost, energy efficiency and user acceptance 
appear to be the most reliable measures of value. Higher lift at 
cruise speeds permits smaller, lighter cheaper wings If landing 
and take-off speeds and rates of climb can be met. Host of the 
new aircraft still on the sketch pads will have smaller more. slender 
and straighter wings demanding all the increased high lift per- 
formance that can be obtained. The figure on the next page shows 
a plan view comparison of the potential dramatic reduction inuring 
size for a 30 passenger transport with advanced high lift. 


I 


I 


Some pay-offs predicted for nev; aircraft applications of the 
best of today’s state of the high lift art are: 

High Speed Single Engine - $ 0 % increases in cruise speed 

and range if C Lniax = 3.0 • 


Short Haul Transport 
(see figure below) 


18 % lower price 

21$ less operating cost 

31# less fuel 


Energy Efficient 
Transport 


18 # less drag (Cr» OY - 3»0, could 
10# less fuel 

use £.0) 


While very large improvements are obtainable with today’s 
best high lift technology, much more can be expected. The tradition 
al two dimensional approach to high lift development should soon 
give way to full three-dimensional viscous flow high lift invest- 
igations made possible by the massive computer analysis and 
corresponding experimental data processing capabilities now in 
existence and being developed. Complete 3-D high lift optimization 
will, for example, result in a smooth spanwise blending of wing 
high lift devices, ( flaps, modified leading edges, etc.) in contrast 
with today's all too typical usage of highly optimized two-dimens- 
ional airfoils chopped into extremely inefficient short spanwise 
segments having large lossy gaps between. 

In summary, the rate of progress in high lift technology is 
high with improvements of £0# range over current operational use, 
economically attractive new families of aircraft have been ident- 
ified in several thorough preliminary design study programs all 
of which show large benefits from improved high lift. 

‘Two areas of possibly useful activity would seem to be: 

1. Safe reductions of the presently required stall speed 
margins may be feasible through the development of 
better stalling behaviour. 

2. "All Lifting" wing and tail combinations will require 
serious development because of the large potential 
gains in high lift. 
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Computer Modeling Nucleotide Base and Base Pair 
Interactions with Amino Acids and Peptides 

John S. Garavelli 
Plant Sciences Department 
Texas A&M University 


This investigation was undertaken at the Extraterrestrial 
Research Division, NASA Ames Research Center, as a part of 
ongoing research into the origin and distribution of life and 
life-related chemistry in the universe. To achieve an 
explanation of the origin of the genetic apparatus of terrestrial 
life it will be necessary to understand the geometries and 
energetics of the interactions between nucleic acids and 
proteins. Nucleic acids must rely on protein enzymes to 
replicate, and in order to synthesize proteins, nucleic acids 
must be transcribed and translated by molecular complexes 
consisting of both nucleic acid and protein components. The 
origin of such an intimately interconnected genetic apparatus may 
be explained through the study of the evolution of hypercyclic 
systems 1 but such an explanation will require a more complete 
description of the chemistry of the interacting components than 
is currently available. 

An additional reason for undertaking this study is that 
recent advances in biotechnology have been absolutely dependent 
on the still fragmentary knowledge of the genetic apparatus and a 
more detailed understanding of the structure and mechanisms of 
this apparatus will be crucial for further advances. Some of the 
processes which have begun to be commercially exploited and which 
involve the interaction of nucleic acids and proteins are DNA 
replication, DNA packaging, RNA transcription, RNA processing, 
tRNA charging, and protein synthesis. 

This research project consists entirely of computer model 
simulations of interacting nucleic acid and protein components 
which are created, manipulated and analyzed with AIMS, the Ames 
Interactive Molecular Modelling System. The work has progressed 
in five phases. The first phase consisted of correcting and 
reanalyzing the molecular structure data base. The available 
libraries of molecular models included nucleotide bases and amino 
acid side chains but these models contained systematic errors 
which were detected when preliminary energy calculations produced 
inconsistent results. The molecular geometries of these models 
were corrected to conform with the work of Arnott and 
Hukins 2/3 and to maintain a uniform set of hydrogen atom bonding 
distances 4 . partial charge distributions for molecules with 
corrected geometries were reassigned with a program which 
calculated molecular orbitals using the MO-CNDO (Molecular 
Orbital by Complete Neglect of Differential Overlap) approach 5 , 
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The second phase of the project consisted of developing a 
procedure for constructing hydrogen bonded molecular assemblies 
from the corrected molecular models. The procedure which was 
adopted employs the addition of phantom atoms at optimal hydrogen 
bonding distances and directions from potential hydrogen donating 
groups. One, two or three hydrogen bonds can then be formed 
between two molecules by superimposing the desired phantom atoms 
of one molecule on hydrogen accepting atoms of another molecule 
using a series of routines available in the MOLECULE program of 
AIMS. Two, one or no degrees of freedom, respectively, remain to 
be determined either by assignment or by minimization of 
non-bonded interaction energies. 

In the third phase a library of hydrogen bonded molecular 
assemblies was produced. The molecular assemblies consist of 
either amino acid side chains or short segments of polypeptide 
backbone which are hydrogen bonded to either single nucleotide 
bases or base pairs already hydrogen bonded in the Watson-Crick 
scheme 6 # 

The intermolecular energy of each of these molecular 
assemblies was calculated in the fourth phase of the project. 
This was done with a routine available in MOLECULE which 
calculated the separate energy contributions of electrostatic, 
non-bonded, hydrogen bonded and torsional interactions 
empirically derived by Scheraga 7 . At this point it became 
evident that some of the assumptions about the various energy 
contributions might be unrealistic and that more work would be 
necessary to delineate clearly which potential energy terms must 
be improved and which terms might be neglected. 

A fifth phase of the project involved initiation of the 
development of a routine to display graphically the results of 
intra- and intermolecular energy calculations. This routine 
which produces a plot similar to those introduced by 
Ramachandran a will be incorporated into the MOLECULE program. 
The table of energy values produced by this routine can also be 
used in information theory calculations of the apparent 
restricted conformation range of flexible molecules and molecular 
assemblies 9 . 
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The Pacer Ins truidentat ions system? when completed? will display 
precise values of mo eh numbers' altitude? and airspeed to the pilot 
a n cl / u r f 1 i H r 1 1 t e « t e 1 1 g i n & o r w h 1 1 e p a c i n & a n a i r c r a f t u n d e r si 0 i n b 
f 1 i g i 1 1 i e 5 t i ri b . P r 0 s e n t r e « u i r e tti e n t s call f o r' t h e firs t u r s i t to be 
installed in an F~1(H airplane* The system will give the operator 
the ability to freeze the display for close observation and to 
select corrected or uncorrected display data* 

The hardware will consist of a two board computer system* One 
of the boards presently under design at NAG A holds the 
mi c roe roccesso r (CPU)? ari thmetic processor ( API) ) ■> volatile and non 
volatile memory? and serial and parallel input /out put (I/O) ports* 

T h c r e in a i n i n b b o a r d in t h e P a c e r s y s t e m will .i n t e r f a c; e w i t h two 
Garrett pressure transducers and the display unit* Provisions will 
also be made to interface with a flight data recorder and/or a 
printer* The display under present consideration is a lieu id 
crystal display consisting of two lines of forty characters each* 

Throughout the entire design? the broader applications of the 
hardware and software have been of high priority* The main system 
board contains the power and flexibility necessary in almost all the 
on board computer applications* The display unit promises to be a 
vers flexible means of presenting .information in a host of projects* 
F 0 r t h e e r <? a s 0 n s ? the so f t w a r e w r i 1 1 e n 1 0 dr i v e tire d :i s p 1 ay? I/O 
ports? and the arithmetic processor has been written as a sot of 
upper level commands to be used by an assembly language programmer * 

To date the hardware arid software have been developing 
simultaneously* The main thrust of this fellowship has been the 
development of the software necessary to utilize the arithmetic 
processor and drive the display* The result of this work has been 
the i-. 1 re at ion of a set of macros and subroutines that may be called 
from an assembly 1 an bus.de program to perform math functions similar 
t o t h o s o f 0 1 .1 n cl on a s o i. e n t i f i c c a 1 c u later* A m a n u a 1 e n t i. 1 1 a d ' A P U 
(8231) FUNCTIONS MANUAL ? 13-Aug-82 A*W* Headley' has been prepared 
ci a c u 111 & 1 \ t i 1 'i & 1 1 1 o s o f t w a r e w i * :i. 1 1 0 n s p e c i f ice 1 i. y f o r t he A P U and t h e 
C P U being us e d o n t h e m < s i 1 1 s w s t e m b o a r d * ( t h c? t a b 1 e o f c* < 1 n t e n t s o f 

this manual is included in the appendix*) The programs can be used 
with no further modifications using the 8085 cross assembler on the 
PDF' 1. .1/34* Included also in this document is an example 
applications program ( PACERo ♦ SRC ) ? which is the program written to 
perform the calculations required by the pacer project* 

The pacer program starts with two 20 bit input* data words from 
the pressure transducers representing the measured total and static 
pressures* The calibrations are then applied to the inputs using a 
polynomial solution* This solution is accomplished by calling the 
P 0 L. Y * f u n c t i o n w r i 1 1 e n for t h e A P U * N e x t t he- press u re r a t i o .is 
calculated using the DIM* function to perform the division* As the 
program continues the mech number is corrected using the table 
lookup and interpolation function (LOKUP*)* This function looks up 
a specific independent variable value in a table and returns to the 
user the correc t i n te r p o 1 a t e d val 1 . is of the dependent variable* The 



calculation'* continue in this Mariner until the values of roach 
number? altitude? and airspeed are calculated* Those values are 
then converted to A8CJX format using the CONMER function and then 
displayed ♦ 

As described above the PA0ER6.SR0 proa rain makes extensive use 
of the macros and subroutines written for the APU and CPU? 

These functions can easily be used to solve necessary equations 
required by future projects using the APU and CPU thus giving this 
wor k i?i a n y f u t u r e applies t x o n s ♦ 


Each function in the manual 
format} 


is described in the following 


name: 


LOKUP ♦ * DPI »QP2»0P3»0P4»0PS»0P6»0P7 
SYNOPSIS * 

The table lookup function will scorch a given table to 
find the closest values to the given independent variable* 
L i n ear intern o 1 a t i o n i s t h e n u s e d t o f i n d t h e d o a :i red 
result ♦ 


DESCRIPTION! 

The operands are defined as follows! 


OPi 

OPS! 

OP 3 

OP 4 

OPS 

0P6 

OP 7 


address of independent variable to be found in table 
address of first independent variable value 
in table 

a d dress o f f i r s t d e p e n d e n t v a r :i a to 1 c v a 1 u e 
in table 

the number of values of the independent variables 
provided in the table 

the address of a scratch memory area <16 bytes) 
a d d r e s s w h e r e t h e p r o gram w .i. 1 1 J u in p if the n u m her 
stored at address OPI is outside of the values 
p r o v i d e d 1 n the table* 

the address where the result wi,13 be stared 


Before using the LOKUP* function the user must establish 
the table of values using the following format* 


OP 4! 

EQU 

V V 1-1 
A A 1 l 



0 P 2 * 

13 W 

1 

kmkkH ? yyyyl-l 

( s e e 

ENTER * ) 

OF' 3 * 

1 

DW 

k x H ? yyy wH 

(see 

ENTER* ) 


RESULTS! 


NOTE* The above groups of data need 
to be in ascending order* 


The results ere stored at address OP?. 
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APU (8132) FUNCTIONS MANUAL 
TABLE OF CONTENTS 


I* I NTRQDUCTION,,,,, ,♦♦»♦**,*♦»»»»»♦,,,♦,♦♦» 2 

II* DESCRIPTION OF HIGH LEVEL FUNCTIONS , ♦ / 


IV, 


V* 


1, 

ADDITION .,,, 


(ADD,) 

» , ♦ , 5 

2, 

APU TO ASCII 

CONVERSION 

(CONVER) 

, ♦ , ♦ 6 

3, 

DIVISION 


(DIV, ) 

* » , » / 

4, 

ENTER NUMBER 

XNTO APU 444444 + ftff 4 

(ENTER, ) 

, , > , 8 

5 , 

EXPONENTIAL , 

»44444444t444444 + + 444 

(EXP, ) 

, , ♦ , 9 


6 , C F GO TO (IF,) 

/ * NAT URAL LOG ,*»«♦♦***♦»»»♦*♦♦♦»»*♦ ( LN ♦ ) 

8* TABLE LOOKUP ,»»»»,♦,»,**,«,,»»,,» (LOKUP,) 

9, MOVE ONE NUMBER TO ANOTHER (MOVE,) 

10, MULT I PL I CAT IN , (MUL,) 

11, POLYNOMIAL CALCULATION , , ,, , * (POLY,) 

12, P 0 W E R ( F' W R , ) 

13, REMOVE NUMBER FROM APU (EEHOV,) 

14 , SQUARE ROOT ( SORT ♦ ) 

1 ti , SUB I RAO I ION »»,*♦♦»»♦♦,»♦»♦»,»,♦»» ( SUB , ) 
1.6, BLOCK MEMORY COPY , ( XFER ) 


,♦♦♦10 
,,,,11 
,,,,12 
♦ ♦ ♦ » *. A 
,, , , 1 5 
,,,,16 
♦♦,♦18 
,,,,19 

. , , ,20 

n •• 

j. 

o '*/ 

4 4-4 4 
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DESCRIPTION OF LOW LEVEL FUNCTIONS 
1 , 


tffrt4-f + + 4-t»t + t + ttf4»fr + tt 


SINGLE PRECISION i FIXED POINT 
DOUBLE PRECISION r FIXED POINT 

FI.. DATING P 0 1 N I »^o 

CONVERSION BETWEEN (1 - 3) 27 


♦»»,,**»»,♦♦♦**,»» 2 4 

'•) K? 
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DESCRIPTION OF FILE STRUCTURE 
1 , 




, , ♦ ♦ ♦ 


2£ 


ARB I.NG ♦ SRC ♦*,,♦,»♦»,♦♦»,»,♦,»♦♦♦♦»♦»»*♦♦♦, ,,,,«29 
2* A K D B L , S R C ,,,,,, ,,,,29 
3, ARFLOT.B RC ,,,,,,,,,,,,*,,,,,30 
4 ♦ A R C 0 N V * S R C ,,,,,,,, ,,30 
S ♦ A R B U B 2 , S R C ,,,,,,,,,,,,,,,,,,,,,,,,,,,<,,,,,,,,,31 
6, CONVERSION OF NUMBERS to/from APU Format ,,,,,,,32 
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ARSING, SRC listing 

ARDBL. , SRC listing 

ARFLOT ,SRC listing 

ARCONV , SRC listing 

ARSUB2 , SRC lifting 

EXAMPLE PROGRAM listing 

NUMBER CONVERSION PROGRAM listings 
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There are three major techniques for measuring pilot workloads 
subjective meaurement, objective measurement, and physiological 
measurement. The greatest success has been obtained with subjective 
ratings where pilots are asked to complete assorted rating scales after a 
flight. However, there is no rigorous theoretical framework to explain what 
such ratings mean, except perhaps for mathematical scaling models that 
contain no axioms related to the content area of workload evaluation. Since * 
it is well-known that pilots (and other humans) cannot provide 
introspective reports with complete accuracy, it is essential to validate 
these rating data with objective procedures. So far, researchers have not 
been able to obtain systematic workload effects based upon objective 
measures of behavior. | 

I believe that substantial advances in human factors will come about 
only when theory from basic research is combined with practical problems 
(Kantowitz, 1981). Objective measurement of pilot workload is a sufficiently 
difficult and challenging problem to serve as a test bed for application of | 
theory to a practical problem where most investigators have been unable to 
find good answers. I came to NASA— Ames this summer to (l)put my money where 
my mouth is and try to solve an important practical problem by approaching 
it from a theoretical perspective, and (2) get my hands on a GAT flight 
simulator, a piece of equipment not available to me at Purdue. 

Despite the valiant efforts of my NASA colleague Sandra Hart who had 
been sending in work orders since last Februrary, the GAT simulator was not 
flying when I arrived at the beginning of the summer. Once the GAT was 
flying properly, we discovered that the signal conditioner used to record 
analog data from the BAT needed work. Making the necessary repairs, 
checking calibration, and adding a four — position switch to the GAT yoke 
consumed almost five weeks. I used two of these weeks to write a computer 
program to control and record the secondary reaction— time task I would use 
as an objective measure of workload since the programmer assigned to the 
project had quit leaving the group temporarily without a professional 
programmer. Another week was spent writing a complex data analysis program. 

By now I was almost as good a programmer as I was during my graduate 
student days. To further polish my dormmant programming skills I sat in on 
EE 611, televised at Ames; this was most useful and an unanticipated bonus 
of my summer. I am grateful to the co-ordinators for informing me of this 
televised course. 

At long last the day of reckoning arrived: We tested Our first 

experimental subject. The reaction time data were exactly what the theory 

predicted. Alas, the GAT signal conditioner became unplugged during the \ 

experiment and so these data had to be discarded because it is essential to 
record both primary and secondary task performance. As of this writing, | 
the experiment is in full swing and data collection will be completed the 
day I leave (if all subjects show up as scheduled and the equipment keeps 
on running). Tentative results (see figure) give grounds for cautious ij 

| optimism. It appears possible to provide an objective measure of pilot 
i workload. However, firm conclusions need await the complete data analyses |j 


:q be performed at Purdue. Additional work, such as expanding the GAT 
-light scenarios and analyzing the micro-structure and instantaneous 
workload, could be done next summer. 

Experimental research with human subjects is always -frustrating, and 
.n this regard despite its superior equipment and technical support 
facilities, AmtiSs is very much like a university. Things always take longer 
:han you plan. Although the chance to use special equipment was what 
wrought me to Ames, there were many other benefits gained. Discussions with 
Bandra Hart and members of her group taught me a great deal about aviation 
and the problems pilots face. Professors Dan Weintraub and Earl Weiner were 
✓isting Ames and we had many helpful conversations about my project. The 
Wednesday evening lectures were uniformly of high quality; this was the 
first time in years that I have paid attention to events outside my own 
discipline and the effort was well worth it. When I return to Purdue I 
shall try harder to attend lectures outside psychology and industrial 
engineering. All in all, the summer was challenging and exciting. It 
allowed me to discover a new area with interesting problems that are 
related to my own theoretical interests in human attention and information 
processing. 

Reference 

Kantowitz, B.H. Interfacing engineering psychology and human information 
processing. In W. Howell t< E. Fleishman (Eds.), Human performance and 
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Implicit Model Following as applied 
to the Forward Swept Wing. 

ORIGINAL 

David Kozel qjt POOR 

Purdue University Calumet 
Assistant Professor of Engineering 

The purpose of the project I worked on this summer was 
to repeat the work done by Grumman on the control law design 
of the Forward Swept Wing. 

The Forward Swept Wing is a control configured vehicle, 
which has some unstable modes. For this reason, a proper 
and robust controller is essential. Grumman used the method 
of discrete implicit model following to produce their gains 
for the control laws. Also, they added an integrator In the 
forward loop to make the system type one. An algorithm was 
then used on the gains so that the output of the forward loop 
integrator, an unobservable state, was not fed back. 

The first step of the project was to develope the 
command augmented state equations (in general form) for, 
the model to be followed and the aircraft. This was done 
for the longitudinal up and away mode. The model state 
equations were derived from a transfer function of. |he model 
taken from a technical report published by Grumman! 

The aircraft state equations were derived partially from 
the output of a linearization program developed by Joseph 
Gera of NASA, Dryden, and partially from augmenting on 
terms representing actuators and a forward loop integrator. 
Afterwards, both sets of equations were augmented by an 
equation representing the form of the pilot input. Since 
the controller was for a digital fly by wire system, the 
equations were then discretized. 

The second step of the project was to develop an 
algorithm to solve the discrete implicit model following 
problem. This was done by transforming the discrete 
implicit model following problem to the discrete linear 
regulator problem without cross product terms in the 
performance index; which was done by using two transforma- 
tions. Then, the steady state solution to the discrete 
linear regulator problem was obtained using an algorithm 
developed by Vaughan". The result was improved using the 
oracles routines. Then, transforming back the results, the 
discrete control law gains were obtained. 

The third step of the project was to develop an 
algorithm that would take the gains resulting from the 
solution of the discrete implicit model following problem. 

Gains that included the output of the forward loop integrator. 

To gains that did not include the output of the forward 

loop integrator, but did include the values of the other states 
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at that time instant and the previous one. A similiar 
algorithm for the continuous case was developed by Anderson 
and Moore. * 

The fourth step of the project was to write a program 
that give'n the augmented state equations for the model 
and aircraft, the relationship between the two sets of states, 
and the desired performance index, it would descretize the 
system, solve the discrete implicit model following problem, 
and transform the resulting gains as mentioned above. 

The following still needs to be completed: Derivation 
of the augmented model and airgraft state equations and 
their resulting control laws, for the other modes of the 
aircraft; Investigation of the capabilities and limitations 
of discrete implicit model following as a design methodology. 
Both of which can be very much aided by the program developed 
this summer. 
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The project began in conference with Hr. Handel and Ms. JoAnu William*. A 
computer literature search (Mod Urn at 5UNY) by the life acioncaa research 
librarian gave S6 pertinent abstracts. The key words were: rat, hybeidoma, 
monoclonal antibodies, calvaria, osteoblasts , osteoclasts. 

The study resolved into three phases: (a) production of hybridomua to yield 
mouse monoclonal ahtibodies against rat bone cells; (b) separation of! the 
bone cells into subpopulationa; (c) Isolation and characterisation of cellular 
antigens Prow the subpopulations. 


l , Mo non Ion a I Antibodies 


Three clones were sought: (a) against bone cells, (b) against a glycoprotein 
surface antigen whigh binds to both Con-A and Lentil lection, (c) against a 
glycoprotein which binds to Peanut- Lectin and Soybean L'.ctln. The procedure 
followed that ofi Kohler anci Milsteln. (1) 


Throe HALU/c mica w‘co In joct.nl l,n, with l x 10' bon** cUn. The bone calls 
had previous l v \vnvp prepared firms n *r*n ,i h*l 1 rat calvaria by Mn JMnn Williams 
and l r. 3iui.lv Holton, fibre*' mice worn* injected i.n. with the Con-A and Lentil 
glvcoprotcin anu throe ”Uh the Poimifc and Soybean glycoprotein. In the latter 
i."o etvwi, 1.H1 ug in ) '0 u) ,y, \ * injected followed in ton days by 10 u g Most. 


Pour davs nfil.ee the last injections the mice wore sacrificed, filmic seen were 
nav»d and tested tor the antibodies. Antibodies morn present in the three 
experiments. 

The mice splc wns w<’ ro removed; snL'nonytnn suparai nd between l' ran ted slides, 
washed and {uispcndod In scrum-free RPMI medium. Cells Prom a mouse myeloma 
line, SP2/o were also suspend >d in serum -free UP MX . Rqual numbers ofi splono« 
cytos and myCloma cells were mixed and polluted togctlvr at 200g x lu min. The 
supernatant was dia carded. 

One ml of! 50% PliG (polyethylene glycol) was added to the pelleted cells for 
one minute at 3?°C immediately fiollo-ed by the addition of l ml of serum-free 
RPMI for one minute. Another ml wan added for a minute followed by seven 
milliliters. The cells were gently pelleted «ud resuspended in 30 ml oP 
complete medium (RPMP, Petal calf serum, pc n ie ill in - s trnp tomy c in , glutamate, 
pyruvate). This volume was a liquated into 3«96wcll microtiter plates and 
incubated overnight at; 37 °0, 10% 00 2 , humtdiPicd atmosphere. Thu following 
morning, the mddium wan changed to HAT medium (Uypoxanthino-aminoptcrin- 
thymidine). In this medium, fusions op myoloma-splenocyte* will grow but 
not myeloma-myeloma or mplenocy Ue-spl ”nocy tc . 11a IP ofi the volume op the medium 

was removed and replaced by an equal volume ofi HAT on days 1,2, 3, 5, 8, 11. Then 
the medium was changed to UT complete medium. After a week the wails wire 
assayed for antibody production. The positive well contents wore transferred 
to 24 well, l ml plates to which had been added feeder thymocytes. 
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the cell* In the 24-well plates worn selectively diluted in order to isolate 
o*»lhj which could ho cloned Giving the desired antibodies. 

ELISA (enzyme- linked immunosorbent assay) (?) 

We preferred thin technique to RIA as our assay for antibody because it is 
less expensive and leas hazardous, 96 well microtiter plates were coated with 
either cell3 or protein. Cells ware added to wells first treated with poly- 
lysine. Tney were then fixed with 0.37a gluteraldehyde and washed witn 0.1 N 
Glycine. Protein was added at 10 ug/ml in 0.1 NallCO,., pll 9.6 and incubated 
2 hr, 37°C, Suparnnfcant* from the wells of the fused oelln ware added to tae 
coated walls and incubated at room temperature for four hours. Tho plate* were 
washed in 0.1% detergent and goat anUl-uouse antibody conjugated with alkaline- 
nhoaolutage was added overnight at room temperature. In the morning the plates 
were washed a no filled with substrata, PNPP, p-nitroobenvl phosphute(colorleas). 
If the enzyme was present, the substrate chan l ed yellow (p-nifcronhonol). 

Colored walls were «? 'i\sid»red positive indicating that the supernatants contained 
the antibody to or nr 'lain. 


II. Surface Antigens 

bon*, cal'**, l :c 10 7 , were solubilised in 2 ml 0.5% NP-40 (nonidet) for 10 
minutes, 4°C. The tubes were, centrifuged 20 min, 3000 rpm and any insoluble 
material discarded. The contents were dialyzed against PUS (phosphate buffered 
saline) to remove the detergent. The dialysate was passed through a Si ml 
Peanut-agarose (Vmg/l ml agarose) column w ich was washed with PBS. When no 
further protein was eluted ns determined by absorbancy at 2110 nm, the column 
was eluted with 0.2 M gn lactose followed by 0.1 M bombs. Protein was eluted 
and was labelled Galactose-plus glycoprotein. The experiment was repeated 
using Soybean lectin-agarose and the same glycoprotein was obtained. 

The solution eluted with IMIS wash from the above columns was passed through 
1 ml columns of Con-A agarose and Lentil -agarose. Each column was washed 
extensively with PUS and eluted with 0.2 glucose followed by 0.1 borate and 
0.2 sucrose. These elutions were pooled and dialysed and protein concentxation* 
determined. They Were labeled glucose-plus glycoprotein. 

The galactose(^ and glucose ( + ) compounds were passed through Sephadex G-200, 

30 ml columns and eluted with Phs. One peak came through for galactose (,+). 

A major peak came through Cor the glucose («•) compound and a lesser peak which 
was a multiple of the molecular weight for the lacge.r one. 

The proteins were further studied on SDS«PAGE (polyacrylamide gel electro- 
phoresis). Purities were found to be singular and their, molecular weights 
determined which were the same estimated from G-200, 

III, ’ Subpopulafcions of Bone Cells (3) 

A fifteen milliliter solution of 20% fetal calf serum, heat denatured at 
56° 30 min, and PBS was prepared. Bone cells were mixed with peanut lectin 
and layered on this column, A huffy layer formed at the top. Afte r 30 min 
in the cold ^ome cells had sunk to the bottom; most had remained in the huffy 
layer. The cells at the bottom were removed and suspended in 10 ml of 0.2 
galactose in cold for 10 min. They wer*» pelleted Iby spinning at 200g x 10 min. 
These cells were labelled Galactose(f) cells. _ 
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This experiment was repeated using lentil lectin. Again cells sank to the 
bottom of the FCS-PBS column after 30 min in cold. These were ?us»cnded in 
10 ml 0.2 M glucose and centifuged. They were labelled glucose (+) cells. 


We are currently comparing histochemically the three cell populations: 
glucose (+), galactose (+), and glucose-galactose (-}. 


We have used our two glycoproteins to immunize mice and are developing 
hybridomas to yield specific monoclonal antibodies. 


We wish to thank Dr. Adrian Mandel and Ms JoAnn Williams for therr 
cooperation and helps and the administrators of the NASA-ASEE Stanford 
University program for permitting me to be a part of this project. 
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In NASA's space motion sickness research meaningful and valid scientific 
measurements must be taken of the inner ear's vestibular system. It is impera- 
tive that external noise or vibrations on the test equipment be kept to a mini- 
mum. Scientific study done on various specimens’ vestibular systems concludes 
that external vibrations in the vertical direction must be kept below 500y^g's 
(5x10~4 times gravitational acceleration) between frequencies of zero to 100 
cycles per second. Preliminary study on a linear accelerator track apparatus 
indicated that an unconventional bearing system would be necessary to keep 
vibrations below the 500y*g level. Therefore, the suitability of air bearings 
became a necessary research area. 

Air bearings fall into two primary classifications: fixed orifice and 
porous. Vibration measurements were taken on three commercially available 
air bearings and in addition to these a series of measurements was taken on 
a roller bearing supported system built and operated by M.I.T. As expected 
the roller bearings gave vertical accelerations far in excess of the 500/^g 
level. The maximum acceleration had a magnitude of 3.6 mg(3600/*g) at 22Hz. 

Hard bearings proved to be definitely unsuitable for our purposes but inter- 
esting for comparison purposes. 

The first of the air bearing systems tested was one in use here at NASA 
Ames. A 3-axis coordinate measuring machine built by Brown & Sharpe Mfg. Co. 
uses air bearings to hover between 0.0003 and 0.0006 inches above a granite 
•slab. This particular bearing had both single and multiple discrete orifice 
arrangements with a system of grooves which distributed the air flow over th® 
bearing interface. A plot of acceleration as a function of frequency was ob- 
tained by doing a Fourier analysis with data from Sundstrand accelerometers. 
Although this bearing produced audible noise, the vibration magnitudes measured 
in the dynamic mode were acceptable as the maximum acceleration was 340 /*g at 
25 Hz. 

A porous type thrust bearing system found on a coordinate measuring machine 
built by Bendix Corp. was the next type tested. As its name implies, a porous 
bearing consists of a perforated surface resulting from a large number of tiny 
orifices, typically 0.1 mm in diameter. This bearing produced no audible noise 
and gave vibration levels well below that of the Brown & Sharpe version. The 
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peak acceleration was 125/^ at 20 Hz. An interesting note to the tests 
run on this system is that some damping occurred in the dynamic mode as 
the accelerations it produced were smaller in magnitude than in the static 
mode. 


The final system tested proved to be the most impressive and promising. 
A custom air bearing shop. Fox International Inc. allowed us to test a 
carriage/monorail setup designed for laser applications. The carriage was 
supported by eight circular discrete orifice bearings: 4 on the sled bottom 
and o two on each side. This particular design had three orifices spaced 
120 apart that fed into a channel system. After two dynamic runs the maxi- 
mum acceleration was 62 /*g at 55 Hz (shown on pg. 3). Therefore, this dis- 
crete orifice, minimum film thickness bearing has a maximum vertical accel- 
eration half that of the porous type tested. This bearing or a similar ver- 
sion of Fox International's appears to be more than satisfactory for our use 
on the Vestibular Research Facility linear accelerator. 
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My work this summer involved the study of absolute stability regions of 
linear multistop formulas in the solution of tho ordinary dlfforontial equation 

dy/dx « f[x,y3. £13 

My specific aim was to search for explicit 1 inoar multistop formulas having 
absolute stability regions that included larger portions of the imaginary axis 
about zero than those regions of previously known oxplicit formulas. I conducted 
the search within the order-classes of formulas for which the order is tho same as 
tho number of backpoints. Starting at the lower orders, I completed the 
investigation through order three and have the program prepared for tho fourth 
order search. The Euler explicit formula Is the only member of tho class of 
order one. 1 will discuss the results for orders two and three. 

The second order class of such explicit formulas have only one free 
parameter. If we choose that free parameter to be tho coefficient of tho 
backmost y-value [call it all, then these formulas can be represented as 
foil ov/s when applied to equation [13. [subscripts are given In parentheses3 

y( n+1 > = [ 1 -al 3 *y { n > + al*y<n~l> + Jc».5*h*C [3+al 3 *f ( n ) + t-l+aI3*f < n-1 )3 £23 

The values of al which allow stability at are those for which -I < al < 1. 

The stability regions collapse from the circle z » e‘° - 1 when al « -1, to the 
line segment on the imaginary axis between i and -i, when al ° 1. This is 
sketched in the sequence of illustrations shown below in the hV*p!ane. 


.0 
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■ Tho third order class of formulas has two freo parameters. When these two 
parameters are taken to bo tho coefficients of the two backmost y-valuos tsay 
al and a23 , then this class of formulas can bo represented as follows. 

y< n + 1 ) » C 1 -a l-a23 *y( n ) + ai*y(n-l) +a2*y< n-2 > + Cl/123*h* 

CC23+5 ft al+4"a23*f(n) + C-16 + 8’‘al + 16"a23‘ , 'f < n-1 > ♦ CS-al+4*a23 *f < n-2 >3 C33 

Tho values of al and a2 which provide stability at h * 0 aro those which 
satisfy the following system of linear Inequalities. 

1- al > 0 
1 - a2 > 0 

1 + al + 2"a2 > 0 r'j 

The stability regions vary throughout this feasible region, however. In regard 
to the measurement of stability on tho Imaginary axis, they are little better 
than ,the second order methods discussed previously. In fact the familiar 
third order Adams-Bashforth formula Cal => a2 = 03 is very competitive. Typical 
stability regions for members of this class are given below. 
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As mentioned earlier, the programming Is ready for tho fourth order 
Investigation, however, time dtd not permit its completion. A couple of other 
Interesting points Wore noted _ dur lug tho summer's v/ork. I will state the first 
as a conjecture. 


Conjecture: Lot L 

y(n+l) = 


be an cxpl 1c 1t_ 1 Incar k-stop, 
y(n-k-l) + b(i)*ftn-1 >, 


kth-ordcr formula of tho form 


i-c 

where k Is even. Then tho absolute stability region of l Is an Interval about 
zero on the Imaginary axis. 


Numerical results for k « 1, 2, 3, and 4 have been consistent with this con- 
jecture. In addition, a sketch of a proof has also boon supportive. I did 
not take the time to pursue a careful proof, although the direction to pursue 
Is evident. This type stability region exists for a third ‘order formula also 
as wo noted In the sketches following the previous paragraph. I did not 
attempt a generalization for odd-orqored formulas. 

The second point of Interest follows from a consideration of the absolute 
stability region for the scheme defined by the pair of equations bel ow. ** * 


Y(n+1> = Cl + zl*y<n> 


y( n+ 1 ) » Cl + g*z3 *Y( n+1 > C53 

In these equations z Is the elgenvalue-stepslze product, and g Is an arbitrary 
real parameter. The stability polynomial for the scheme above Is 


pCg,zl = 1 + Cl+g3*z + g*z**2. 


C63 


The following result Is established. 

* » * » 

Theorem: The closure of the absolute stability region of the scheme C53 

contains the boundary of the absolute stability region of Euler’s explicit 
formula for 3 < g < 1. 

Proof; The absolute stability boundary of Euler's explicit formula is 
z = e; - 1-. Consider pCg,e t6> - 13 , where p is taken from C63 . The magnitude of 
. pCg,e'°- 13 Is the magnitude of’ the quantity * * ;• 


g*e* tf + 1 - g. 


C73* 


Of interest are those values of.g for which C73 Is contained in the closed unit 
disk. Since the unit disk, g*e‘° , and C73 are each circular with center on the 
real axis, we may reduce our concern to 0=0, and ir . When & <= 0, C73 
reduces to 1 for all g. When 0 a 1T , C73 becomes 1 - 2*g which is in the closed 
unit disk for 0 < g < 1. 


*** From oral communication with Harvard Lomax. 
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NONLINEAR SOLUTION TECHNIQUES FOR THE DISCRETIZED 
STEADY-STATE EULER EQUATIONS 

Charles W.Schelin 

The Euler equations amount to a system of conservation 
laws that may be used to describe inviscid fluid flow. Many 
numerical schemes have been developed to solve these equations 
using time step methods. If one Is concerned only with the 
steady state solution the iterative steps toward this solution 
need not be time accurate. Thus it may be possible to obtain 
the steady-state solution more rapidly using nonlinear schemes 
whose iterative steps are in some sense independent of time. 

Dennis Jespersen of NASA-Ames has successfully applied 
Newton's method to the steady-state Euler equations using a 
Gauss-Seidel multigrid scheme to solve the linear system for 
each Newton step. I began this summer's work trying to speed 
up his scheme. 

After a few week's work, I was able to reduce the comput- 
ing time required on some of the intermediate Newton steps by 
Iterating only on those nodes where a sufficiently large change 
was indicated. Regular Newton steps were still necessary, how- 
ever, at the beginning and to finally obtain a converged solution 

In the hope that a faster scheme could be developed - I 
spent the remainder of the summer analyzing the nonlinear dif- 
ference equations that arise when the steady state Euler equ- 
ations are discretized using the flux-vector splitting of 
Steger and Warming [1] and the second order one-sided differ- 
ences used by Jespersen. These equations are linear combinations 
of the generalized flux-split vectors E and F evaluated at grid 
point (i,j) and up to eight neighboring points. Letting d,u,v,c 
represent density, horizontal and vertical components of velocity 
and local sound speed, respectively, one finds that the vectors 

E = (el,e2,e3,e4), and F=(f l,f 2,f3,f4), 
satisfy 

el = d*(.8*wl 4- w2) 

e2 = u*el + d*c*w3 

-\ e 3 = v*el 

e4 - u*e2 + v*e3 - .5*el*(u**2 + v**2) + 2.5*w2*d*c**2 

fl » d*(.8*zl + z2) 

f2 - u*fl 

f3 » v*fl + d*c*z3 

f4 » u*f 2 + v*f3 -.5*f l*(u**2 + v**2) + 2.5*z2*d*c**2 
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where 

wl « (u +/-abs(u))/2 

w2 « u +/-(abs(u+c) + abs(u*c))/2 

w3 “ c +/-(abs(u‘fe) - abs(u-c))/2 

and zl,z2,z3 are obtained from wl,w2,w3 by replacing u with v. 
Using these properties, I considered a point Gauss-Seidel solution 
scheme. 


♦ 

A point Gauss-Seidel scheme amounts to solving for the var- 
iables djU,v,c at point (i,j) using the latest values of these 
variables at neighboring points. Such a scheme thus requires, sol- 
ving a nonlinear system of the form 

gl(d,u,v,c) « bl v 

g2(d,u,v,c) « b2 ,! 

g3(d,u.v,c) « b3 

g4(d,u,v,c) » b4 

at each grid point. Taking T = d*c, if bl > 0 and b4*bl > .5*(b2**2 
+ b3**2) , I have shown that this system has a solution. When u .< c 
one finds 

v - b3/(2*X'*T + bl) 
u = b2/(2*Y'*T + bl) 

c - .8*(Y'*abs(u) + X'*abs(v))*T/(bl - 2*(X'+Y')) 

i 

* * * ♦ ♦♦ 

where X' and Y' are derivatives of the coordinate transformation 
and T must satisfy g(T) = 0, where 

. • g(T) « 5*(X'+Y')*T*c**2 r b4 + b3*v + b2*u - .5*(b2**2 + b3**2) 

In 1 the supersonic case the formulas for u.and c are more •. 
complicated, specifically 

u = (bl*b2**2-( 1.6*X'*Y'*abs(b3)+4*b2*X'**2)*T**2)/h, 

where 

h » (2*X'*T+bl)*((5.6*Y # *T)**2-2*X'*bl*T+bl**2) , 

and 

c « (b2-bl*u)/(2*Y'*T) . 

In this case, g(T) is altered by replacing 5*(X'+Y')*T*c**2 with 
5*(X'*c+Y'*u)*T*c. If T' is selected so that u = c, and if g(T') < 0 
then the system admits a unique solution T in the interval (T',R) 
where u < c and R =» bl/(2*(X'+Y'))- If g(T') > 0, the system ha© a 
solution in the interval (0,T') where u > c, but none with u > c. 
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I have written a FORTRAN code for this scheme which I am 
taking back with me to UW-LaCrosse for further refinement. Thanks 
to Dennis Jespersen and the other fine people at NASA-Ames this 
summer's experience has opened an exciting area of research to 
me. With continued work I hope to obtain results that are worth- 
while computationally. 
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This project concerned the computer modelling of a viscous incompressible 
fluid flowing past a circular cylinder. The flow is assumed to have constant 
velocity at infinity. The problem reduces to that of a two-dimensional flow 
past a disk. 

Properties of the flow are determined by the Reynolds’ number, which is 
defined to bo „ Ll w * H 

K - — ' 

where is the speed at infinity, j£ is the diameter of the cylinder, and 
7 J is the viscosity. The velocity vector of the flow satisfies the Navier- 
Stolces equations: 

-t- u- Vu - V V\ - Vp 

V” ■ a — c 

where p is the pressure. A computer program to solve these equations numerically 
has been developed by Unmeel B. Mehta of NASA-Ames [1]. A listing of his program 
may be found in [43* 

There is a critical Reynolds' number Rq such that for R< Rq , the 
flow eventually converges to a steady-state flow. The flew is stable, meaning 
that any small perturbation quickly dies away. When R'> R^ , the flow becomes 
unsteady, with eddies constantly forming and being carried downstream. It has 
been experimentally determined [3] that Rq = 48 + 1. The generally accepted 
explanation for this is through a bifurcation of the solution td the Navier- 
Stokes equations. There is, presumably, a steady-state solution for R */ R Q 
but it is unstable. 


My project, suggested by Hurray Tobak of NASA-Anes, was to obtain some 


numericajl evidence for this explanation through computer simulations. The 
Nnvier-S tokos equations, written in terms of the vorticity , are: 

| L 7 ? T -- v V'V + Vf y- V 7 'l|- 

The linearised equations arc: 

f { v 2 \ -< W'-f ' V}*y\ t v-s *vy 

This leads to the eigenvalue problem: ^ ^ " "0 \7 ^ i’ X V* ^ 'V’\ 

Using values for l|y obtained from Unmeel Mehta’s program, it is possible to 
solve for the eigenvalues }\ on a computer. If R < 48, the eigenvalues should 
have negative real part. If R is close to 48, there should be a pair of 
eigenvalues 

^ = a -H b i , (a,b real) 


Where a is close to zero. Then b will be related to Che experimentally 
observed frequency (see [3]) at which the unsteady solution oscillates. 

I used a 128 x 84 grid on an annulus, so the matrix to be diagonalized is 
10584 by 10584. I have prepared a computer program to do this using the 
following strategy: By reducing to a coarser grid, say 32 x 24 , it is 

practical to use the QR algorithm to compute all of the eigenvalues and 
corresponding eigenvectors. The only eigenvalues of interest are the ones with 
real parts near zero. The accuracy of these can be improved by going back to 
a finer grid and using an iterative method such as the Wiclandt method (see [2]). 
Preliminary results indicate that this strategy is feasible, but so fax- 
numerical instabilities have prevented getting results accurate enough to 
compare with experiment. I hope to be able to complete the project next summer. 
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I have made considerable progress in the effort to model magnetic braking 
in protostellar clouds. I have adapted the 3D numerical hydrodynamic computer 
code of Tohline (1978) to include the effects of magnetic fields. 

I had earlier modelled magnetic clouds in 2D using a vector potential 
representation of the magnetic field. This commonly used method is not appli- 
cable to non-axi symmetric field geometries, so that another approach had to be 
devised. This method, called the "FLUX" method, relies on the principle of 
conservation of magnetic flux through any surface comoving with the conducting 
fluid. Individual fluxes through elemental areas defined in each computational 
cell are used to represent the field. The advantage of FLUX is that it is very 
similar to the "donor-cell" method of the Tohline code for transport of the 
various physical quantities such as mass density, momentum density, etc., and 
thus is well adapted for inclusion in this code. 

FLUX had earlier been tested in a 2D version and proved satisfactory. This 
summer FLUX was used in a 3D version, which required much more complicated coding, 
thus requiring extensive "debugging." The aim was to reproduce an earlier 2D 
simulation' of a collapsing magnetic cloud. This provides & crucial test which 
the 3D code must pass in order for any later results to be trusted. 

The accompanying graph shows the results of this test run. Shown is the 
run of central magnetic field strength versus central gas density during the 
collapse of an initially uniform, spherical gas cloud. The evolution follows 14 
an irregular early course followed by a steady collapse with “central^ ^central 

(This behavior is described in detail in Scott and Black (1980)). The solid 

line shows the earlier data derived from the 2D simulation. The "X's" are data 
points from current 3D simulation, also starting from an initially uniform 
spherical cloud. As can be seen, the 3D data are almost identical to the 2D 

data. Unlike the 2D code, the 3D one is not constrained to remain axi symmetric. 

However, over the whole evolution, no spurious non-axi symmetry developed to 
better than one part in 10 IU using a measure that is essentially ( /p ) 
non-axi symmetric. Because we are interested in precisely the detailed rotation 
history of the clouds this excellent performance on the test case is very 
encouraging. 

Having successfully reproduced earlier 2D simulations, the 3D code is now 
ready to be used to attack problems with both rotation and magnetic fields, in 
particular magnetic braking of the protostellar clouds. 
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During the course of the program, two topics related to inviscid 
transonic flow were investigated. The first involved determining the 
possible existence of limit cycles, in unsteady potential solutions. and 
the existence of multiple steady Euler solutions for two dimensional . 
transonic flow over airfoils. The second topic involved the inclusion 
• of convectihg vortices in three dimensional transonic potential flow. 

In the first investigation, computer runs were first made with an 
unsteady small disturbance. computer code, LTRAN2, developed by Peter 
Goorjian and William Ball haus of NASA-Ames(^ These runs were made to in- 
vestigate the time-dependent behaviour of solutions for flow similar to ... 

* ..(21 

that in which multiple steady solutions were found. 1 1 The steady solutions, 

' for symmetric airfoils at zero angle of attack, consisted of a symmetric - 
zero lift one and two mirror-image non -symmetric lifting ones. The first 
study showed that the symmetric solution was the unstable one. Fig.l . 
shows how lift increases with non-dimensional time after a small transient, 
for one quarter time period 1 ', perturbed the initial, symmetric flow. This 
result has implications for transonic buffetting where the physical flow 
does not reach a steady state but continues to oscillate. Our results in- 
dicate that if due to boundary layer or other physical effects, the non- 
symetric solutions are unstable, then the entire flow will be unstable. 

The time units in the figure, correspond to one period (observed experi- 
mentally) for buffetting for a similar airfoil. The slow growth of the 
computed instability compared, to the physical phenomenon indicates that 
another physical effect not included in the calculation, such as boundary 


i 
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layer growth, may be responsible for the observed rapid oscillations. 

Also, self - sustained oscillations were not seen in the computed re- 
sults and the unsymmetric solutions appeared ‘to be stable, futher in- 
dicating importance of other effects such as boundary layer growth in 
buffetting. 

A second part of this study involved the steady Euler equations. ■> 

A series of computer runs was made using a two-dimensional Euler Equation 
solver developed by Thomas Pulliam and Joseph Stegeri 3 ^to determine 
whether non-unique solutions exist for the Euler equations. A systematic 
method for investigating" this question involves determining the amount of 
lift induced by a small angle of attack as a function of Mach number. 

Results inRef. 2~Tndicate that this quantity becomes singular at the 
bifurcation point where multiple' solutions appear. Preliminary results 
indicate that, for Euler solutions, this quantity is not singular, at * 
least until Mach numbers are reached for which the code fails to converge. , 

This convegence failure may indicate lift divergence. Further investigation 
of this question Is required. 

The second main topic;' involved embedding a convected vortex in a 

’ ' * (A) 

three-dimensional rapidly converging potential flow code'. The test case 
computed involved convection of a line vortex past a wing, where a weak 
shock was induced. The rapid convergence of the code was not affected by 
the presence of the vortex. 
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ISOTOPIC RATIOS IN SPACE DUST AND ATMOSPHERIC AEROSOLS 


Alien B, Tucker 
Physics Department 
San Jose State University 


At MeV energies, individual nuclei can be counted and the size 
of the pulse produced is a measure of the nuclear charge. Thus mass 
snectrometry with nuclear particle accelerators is capable of very 

—1 2 —T r 

high sensitivity (assays at 10 - 10"* ? ratios of desired-to~ 

background isotopes have been reported) with small samples (about 
10^ atoms are needed in a sample for a measurement with Vfo counting 
statistics). This summer I have done preliminary work on application 
of this new technique to two NASA/Ames research programs. 

First is analysis of "space dust". As it moves about the sun, 
the earth sweeps up interplanetary dust. In contrast to meteors, 
small particles (about 10 microns) can radiate the energy associated 
with their cosmic velocities rapidly enough to survive the capture 
essentially unaltered. Collected in the stratosphere on U-2 air- 
craft with minimum terrestrial contamination, they are a record 
of the composition of the early solar system. I am developing a 
proposal for isotopic analysis of selected light elements at the 
new NSF-U. of Arizona accelerator mass spectrometry facility. A 
preliminary test of ion source output from a 20 mg sample of the 
Murchison meteorite was conducted in August. Beams of 10 parti- 
cles per second of the abundant stable isotopes were produced. 

A second application is measurement of cosmogenic Be~10. Radio- 
active Be-7 and Be-10 are produced in the stratosphere by cosmic 
rays and adsorbed on aerosol particles, With half life of 53 days, 
Be-7 activity can be counted by conventional gamma-ray spectroscopy. 
The 1.6 s 10 6 yr half life of Be-10 makes decay counting impractical. 
However the Be-10 inventory in samples of the size collected by 
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the U-2 aircraft can readily be counted by accelerator mass spectro- 
metry. Since the short-lived Be~7 activity will tend exponent i ally 
toward saturation (decay rate equal to production rate), ratios 
of Be-7/Be-10 would be a measure of stratospheric residence times 
for the air masses sampled. The new NSF-U. of Arizona facility 
has not yet measured Be-10. However labs at IT. of Pennsylvania, 

Or say, and ETH Zurich (with which 1 have contact) have experience 
with this isotope. A proposal for Be-10 counting of aerosol samples 
is being written. 
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Electron Extraction and Attendant Processes 
Induced by Ions Incident on Crystalline Metallic Surfaces 


I. Harold Zimmerman 

Assistant Professor of Physics 
Clarkson College 
Potsdam, New York 13676 


The investigation undertaken this summer is the first stage of a rather 
ambitious program of research aimed mainly at understanding how electron trans- 
fer takes place when a positively charged ion is neutralized during encounter 
with the surface of a metal. Eventually I hope to include nuclear recoil, lat- 
tice dynamics and some electronic structure effects in the model, but for the 
present a very simple approximation to the system is being employed. This still 
should suffice to demonstrate the basic neutralization process and to give some, 
indication of how free electrons can be produced by the ion impact. Understanding 
these two phenomena, electron transfer to the ion and the production of free 
electrons, is the principal goal of the present stage of research. 

The metal is represented by an unidifferentiated, positively charged 
"jell i urn" that occupies the left half-space. The negative x-axis lies with- 
in this jellium, with the origin at the jellium surface and the positive x-axis 
extending into the vacuum. The surface thus coincides with the y-z plane. 

Ion bombardment is represented by a point charge of magnitude +q moving 
toward the surface with a specified velocity along the x-axis. The cloud of 
conduction electrons associated with the metal is modeled in accordance with 
the density functional formal ism* as a charged quantum fluid whose charge 
density is equal to that of the jellium deep in the bulk. Nearer the surface, 
the electron "fluid" decreases in density and extends into the region of positive 
x owing to its self-repulsion. Within this model there is clearly no possibility 
of representing lattice dynamics; the lattice has been smeared out into a uniform 
background of positive charge. This, in turn, means that the trajectory of the 
incident ion cannot be ascertained from the dynamics of the model. It must be 
prescribed. That will do for the present, however; our first interest is in the 
electronic response to the impact. Our simple model lets us address several relevant 
questions: (1) Are collective electron effects important? (2) Is the neutral- 

ization accomplished principally by Auger-assisted transfer to the incident ion's 
neutral ground state, or by quasiresonant transfer to excited levels of the ion? 

(3) If Auger transfer dominates, does it manifest itself by means of convoy 
electrons traveling in train with the receding (neutralized) ion, by means of 
excitation within the electronic fluid that remains associated with the metal, 
or by some combination of these? 


0 


2 


The energetics of (quasi-) resonant and Auger-assisted electron transfer 
are sketched in Fig. 1. (Quasi-) resonant transfer can easily be a one-electron 
process, but Auger transfer must involve at least two electrons and possibly more. 
A reasonable scenario for Auger-assisted transfer involving the production of 
both convoy electrons and internal excitation within the metal is indicated in 
Figure 2. This figure also suggests how the results of these calculations might 
best be displayed: by means of computer-generated motion pictures portraying the 

evolution of the electronic fluid under the stimulus of the ion impact. The 
computer program is therefore designed to facilitate eventual translation of 
the output data into graphic form. 

Initially, a 1-dimensional approximation is made whereby -he complexity 
of the calculation is greatly reduced as is the time needed for a run. These 
reductions let us try out a number of computational schemes so that the most 
suitable approach will be evident before any large investment need be made. 

Within the program the electron density amplitude function is represented 
using the finite element method in conjunction with the Raleigh-Ritz-Galerkin 
variational procedure 2 . The functional is written as 

j = Jcitpx [ af t - h] , 

where is the density amplitude function and H is a local approximation to 
the multi -electron hamiltonian operator^. The ranges of the integration are 
assumed to be such as to justify the assumption that Tjf never changes on the 
boundary during the times of interest. Within the region of interest, 
changes in such a way as to render J stationary. 

The problem is nonlinear because II depends on $. This leads to a number 
of challenging difficulties that have taken up most of the summer's effort. 

These seem now to be largely in hand, and the program can profitably be trans- 
ferred to the Clarkson computer system when I return there. By next summer 
I expect to have finished work on the present 1 -dimensional model system and 
to hsve the full 3-dimensional calculation ready to run. Treatment of nuclear 
dynamics should also be ready to begin at that time, represented by the ad- 
dition of a lattice core ion positioned at the impact site. 
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Dictyontoliuin disco ideum , an Ontogenetic Modal Far the 
Study of Bio-Gravitational Effeata an Cellular to Organ ioraiu 
Programmed Development, Including Senescence nnd Death 

Morris Aaron Bonjoniinson 

Lang Inland University College of Pharmacy & Health 
Associate Professor of Public Health 


Km Dhj ledge of the unique life cycle end o trot eg, in evolutionary position of Dictyo - 
ntelium diEicoideum prompted me to suggest it as a candidate for free fall ex- 
perimentation in a response. I submitted to the first solicitations for life soiencB 
experiments to bo flown aboard the Space Shuttle. 


This organism begins its career, upon emerging From a spare, os an independent, 
actively feeding and dividing myxameba. Its food consists of bacteria. After 
feeding and cell division have simultaneously ceased, certain myxomobae send 
out pulses of chemical signal (hormone) which, upnn contacting other arnobaa, cause 
them to elaborate pulses of the same chemical messenger (cyclic adenosine mono— 
phosphate) arid tn migrate in the direction oP the primary signal. These amstaao 
become secodary sources, their target cells, in turn, become tertiary sources, 
and so on (1). The migrating cells clump together to form an aggregate which under- 
goes a series of morphogenetic movements with concomitant cellular differentiation 
to form an upright structure erected at a 90° angle to the substrate. This stage 
is colled a papilla. The papilla then lies down upon the substrate farming a 
grox (slug) which migrates far a period of time. It then becomes stationary and 
undergoes sub-cellular and grass morphological reorganization to n sorooarp. This 
consists, mainly, af a slender, tapering stalk bearing at its apex a clump of 
spnren. These are disseminated and repeat the life cycle. 


All of these events are under genetic control. Gene action is modulated by in- 
trinsic and extrinsic environmental factors. The substrate (agar, in the labora- 
tory) should be thought of no an integral, part of the organisms' functional anat- 
omy since antagonistic, hormone -like substances reach their target cells by diffusion 
through the interstices of the agar lattice. Depending upon genome and environ- 
ment, the Ufo cycle is completed in from five to seven days. 


My experiments have been done on the ATCC-Raper and CB wild type strains of D. 
disenideum . I have used the standard tools far determining gravity effects ,~t ho 
clinostat end the centrifuge. There is an on-going debate as to the validity 
of considering rotation on a clinostat a means of nullifying tine effect of the 
gravity vector (2, 3, 4, 3,(5, 7). Salisbury ( A), however, based upon observations on 
the space flown and Earth experiments, states, "the clinostat simulates free 
full at the cellular level in many important respects." There is, on the other 
hand, little doubt that increased gravitational force results from centri- 
fugation. The results of my clinostat experiments are, in summary; (1) Rota- 
tion of cultures in a range of RPM's with the substrate parallel to the gravity 
vector, decelerates development and produces a larger number of organisms when 
compared with static horizontal controls, regardless of inoculum stage (spore 
or myxnmuba) and genotype, bearing in mind that LID 1 s rate of development is, 
under standard condition;?, slower than ATCC's. (2) Rotation of cultures with 


the substrate perpendicular to the gravity vector decelerates development arid 
often shows fewer organisms than do horizontal controls* (3) Inverted, static 
plates, with agar surface facing downward, show accelerated development relative 
to static cultures incubate'' with the agar surface facing upward. Inverted and 
rotated plates develop even more rapidly. (4) Non-rotating cultures maintained 
with the agar surface parallel to the gravity vector, usually show decelerated 
development and an increase in numbers of organisms when compared with horizontal 
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.literature included descriptiona a I' BX05ATELLITE IX, DISCO' 
SMYtAB 3,4, GEMINI 3,11, NERU 1, APOLLO -SO YUZ, APOLLO 16,1V, 
HOSMQS repnrto, ate. 
ranging Pram microbes: 


'R XII, XU1 I , 

COSPAR volumeo, 


A apnea flight effect which has bean nerheu in diverse bint a , 
tn insects i;n plant material to fish eggs, in a change 


in the rata nf the genetically prugrammerj life cycle relative to that of Earth 
controls. D. dinnnltimim 's life cycle: emergence, feeding and growth, cell 

differentiation, reproduct ion/ounHSconcB , death, lands itself tn studios of 
developmental rates. Last year, at Ames, in preliminary stud inn with the fluo- 
rescent vital dye, ^-rioetnmidu-V inathincynnouti, Iberia 2,2' dinulfonln acid (SITS), (16) 
I was able tn Unmnnn brute a lack of toxic of feu bn throughout the Dictyasto.lium 

life cycle and to differentiate, boned on the dye's staining properties, between 

living and deed I), dinenidnum cello, BITS provides sensitive, accurate meann far 
monitoring the entire life cycle of the organism. Ilutz (17) lino pointed out the 

analogies between the developmental program in D. dlacoideum and that in de- 

veloping embryos. Thu time span in Dlctynntel him is, in mast canon, shorter. 

Pollard (11), based on bin calculations, believes that studies an developmental 
systems greater than 1 mm would be both promising and interesting. D. diocniduum 
is a developmental system wlriuh achieves this minimum oi y.c and beyond, goes 
through a complete developmental life cycle in 5-7 days which in similar to 
divergent life forma, and includes, on analog of embryogenesis in its programme. 

It is, therefore, an ideal candidate fur experiments in BPACELAB. The instrumental 
parameters for integrating D. dincoideum experiments into the BPACELAB configu- 
ration can be easily determined through use of a device such as the Plant Carry- 
On Container (PLUG) in its designated rule in a test similar to tine HEFLEX Bio- 
engineering Tea i; (HITT) (IB). 
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DISPLAY RIIQU I RKHKNTS PCJR REMOTELY PILOTED VEHICLES 

Barbara A» Bern ah o 
Florida Institute of- Technology 

A fair amount of at:, ton Li on at Dr yd on Flight Research 
F&ci 1 i ty has foe. used on the functional s operational , and handling 
char scfceri ski i:s of remotely piloted vehicle systems (RPV) • In 
such it system, a pilot, controls an aircraft from a ground based 
cockpit in t>i ther a full manual or coupled control mode. 

The information available to the? pilot for flying 
•maneuvers such an spin initiation,, spin recovery, and approach and 
3 and i ops is limited to instruments and a TV mor»i ter which displays 
a vidua image transmitted from a camera mounted on the aircraft, 
this presumably provides a pilot's eye view of the world. 

To illustrate,, one such system is the HI MAT, acronym for 
highly maneuver abl e aircraft: technology. It is a 3/8 scale model, 
flown from a ground based cockpit. Such a system raises a number 
of serious issues with respect to human information processing, 
manual control , and pilot workload demands. On the one hand, 
unlike a simulation there is something at stake. While not iif e 
threat eni ng, a serious error* could crash the aircraft not to 
mention crush the pilot's ego. On the other hand,, flying an RPV is 
something like doing an impoverished simulation, eg. in a fixed 
based simulator , since many of . the cues present in flight are 
unavailable in the ground based cockpit. 


PERCEPTUAL CPUS I DERAT 1 OHS 

i„ Cockpit. Layout — Since there is no si 
the cockpit instrument layout is spread out,, thus inc: 
necessary eye and head movements of the pilot,, 


e constrai nt , 
easing the 


2. Vestibul. ar/Propri ocept.i ve Cues Cues to orientation 
in space are absent. This in itself is not unmanageable since 
visual cues can often override or substitute for these. However , 

:i n RPV flights visual cues are also i mpover i shed . 


!. . i m i t e t;l V i s u a 1. F i e 1. d 


constrained to a 19 " TV monitor presenting a 


i i e p i 1 at ' e v i sua 1. f i e .1 d i s 
?4 degree 

horizontal and .1.4 degree vertical, field of view. This poses two 
kinds of problems. A) Specifically in approach to landing, it can 
distort perception of position with respect to the runway when 
the aircraft is banked or in crosswind condi t 
heading is corrected for crosswi nd the runway 
completely outside of field of view. The b& 
approaches where the horizon may fall outside 
B> A limited field of view also r 
el imi nation of peripheral vision which supplies important motion 
cues that in turn contribute to perception of altitude and 
ai. i "speed » ’I hi © happens because objects appear to move more slowly 
across the visual field than nearer ob iects. The pilot typically 
uses t.Iiis subtle information in his assessment, of altitude. 

Current research at Dr yd on includes a wide angle video 
system provide 1.60 degrees field of view and limited resolution 
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GRAVITATIONAL EFFECTS OF THE ASSOCIATION,NGC 206, ON THE DISK OF THE ANDROMEDA 
GALAXY 

Gene G. Byrd. Associate Professor of Astronomy, University of Alabama 

The galaxy M31 h^s a two-armed spiral pattern in Its disk. Recent 21 cm 
observations by Unwin 1 2 show the hydrogen gas structure of the disk of M31 with 
improved resolution., The main spiral arms can be seen in the observations, but 
there are a number of irregularities. In particular, a steeply pitched "armlet" 
extends along the southwest semi-major axis of M31's tilted disk between the two 
major arras (see figure 1). 

My project has been an attempt to explain this structure. Photographs of 
M31 show a large star cluster, NGC206, at the inner end of the armlet (figure 1). 

My hypothesis was that the armlet and other nearby disturbances are caused by the 
gravitational effect of NGC206 . I intended to quantitatively verify this using * 
computer simulations of NGC206 orbiting in the disk of M31. 

Last summer was spent on two approaches to this project. I worked on developing 
a program for a first, simpler, approach to the problem. The second approach, 
more realistic but more complicated! involved modifying a program by R. H. Miller 
and B. F. Smith (my NASA colleague). During this last winter, I continued work on 
the simpler approach and submitted the first version of an article describing 
the results to The Astrophysical Journal , 

During this summer, I have completed the simple approach computer runs and 
written the final version of the article which has been accepted for publication in 
the Jan. 15, 1983 Astrophysical Journal . The results of the simpler approach showed 
that less extensive modifications than previously planned would be sufficient in 
Miller and Smith’s program. I have finished these modifications and am doing 
computer runs for the more realistic simulations. The article results show that 
NGC206 can create the armlet by its gravitational effect on material in the disk 
of the Andromeda Galaxy. These results have wider implications , supporting 
ideas that gravitation is the dominant mechanism for generating spiral arms in 
galaxies with multiple arm patterns and indicating that mechanisms involving chain 
reactions of supernova explosions and star formation are not important in spiral 
arm formation. 

I will take the results of the more realistic simulations back to Alabama 
for analysis. Discussions with P. M. Cassen, another Ames astronomer, have raised 
the possibility that gravitational armlet forming processes may have been important 
in the much smaller scale process of forming Jovian planets in our solar system. 

I will investigate this interesting possibility further when I return to Alabama. 
Finally, I have atAlabama for furter analysis the results of simulations of- 
collisions of elliptical galaxies by Miller and Smith. 

1. Unwin, S.C. 1980a, Monthly Notices of the Royal Astronomical Society vol 190, p. 5511 

1980b, Monthly Notices of the Royal Astronomical Society vol 192, p.24;j 

Baldwin, J.E. 1981, The Structure and Evolution of Normal Galaxies (New York: t 

Cambridge University Press) S.M. Fall and D. Lynden-Bell eds. p.145. 

2. Miller, R. H. and Smith, B.F. 1980, Astrophysical Journal vol. 235, p. 421. jj 
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Determination of Asslmilatory Quotients 
for Chlorel la pyrenoldosa 71 1-05 
Grown in a Closed System 


Arthur P. Carroll 
Associate Professor 
Department of Science 
College of Santa Fe 


In order to sustain extended manned missions in space ( missions over five 
years), it will be necessary to develope regenerating life support systems . The . 
CELSS program ( Closed Ecological Life Support Systems ) at NASA/Ames is involved 
In research to that end. 

One of the most basic problems in sustaining life in closed systems is the 
need to regenerate gasses. Animals or heterotrophs consume oxygen and produce car- 
bon dioxide, while plants, autotrophs, consume carbon dioxide and produce oxygen. 
In a closed system gas balances could be maintained by combining autotroph and 
heterotroph. Unfortunately the ratio of gas consumed to gas produced Is not the 
same for each. The quantitative expression for the CO? consumed and the 0 2 producd 
by plants is known as the Asslmilatory Quotient ( A ) 

A great deal of information regarding the Asiimilatory Quotient of plants, 
especially the green algae, exists in the literature. However, most of these stud- 
ies occurred in open, flow through systems rather than closed systems. In order 
to study the effect of closed systems on the A , NASA designed a closed system 
to accommodate growth of both an heteroph and $n autotroph .('Figure 1 ). 

Before the system could be fully utilized, it was necessary to test its 
reliability and reproducibility. Last summer experiments were performed using 
the green alga Cioreila pyrenoldosa 711 ~~05 . The A of the organism was determined 
in the closed system under conditions of d if ferent^O? starting concentrations. 

The initial data indicated a variation in the A q at different concentrations. 

This summer experiments were performed to verify the results of last summer. 
It was necessary to determine the varience in results obtained at constant condi- 
tions. Six experiments were performed to determine the A at a starting CL con- 
centation of 20%. ( Table 1 ). It quickly became apparent that reproducible 

results were difficult to obtain. The initial experiments were performed at a 
reduced pressure ( 650 mm Hg initial ) in a Helium background. Increases in pres- 
sure would occur during the experiments. These increases, if due to leaks in 
the system^would account for the variations In results. An experiment was perform- 
ed in which the starting pressure of 650 mm Hg was used and the gas In the head 
space was monitored by means of gas chromatography. The presence of Nitogen in 
the closed system would indicate leakage. Nitrogen was found and the quantity 
present agreed with the increase in pressure. Two more experiments were run with 
starting pressures of 760 mm Hg. The A of these experiments showed closer 
agreement. - q 

Because of the existence of leaks the data produced can give erroneous 
results. The system must be used with caution at reduced pressure. The high rela- 
tive humidity in the system can cause the vacuum fittings to leak gas from the 
environment. The system does seem to be reliable when used at atmospheric - pres- 
sure . Experiments done at 7&0 mm Hg, when monitored for the presence of nitrogen 
gas, showed none. 



TABLE Is Assimi latory quotients for Chlorella pyrenofdosa 711-05 grown in a 
closed system at initial oxygen concentrations of .20% and carbon 
dioxide concentration of 5%. 


Experiment # 

V 

Starting Pressure 
( mm Hg ) 

Pressure Change 
( mn Hg ) 

1 

0.72 

V. 

650 

70 

2 

0.81 

650 

64 

3 

0.68 

650 

83 

k 

0.59 

760 

6 

5 

0.57 . 

760 

6 

6 

0.64 

650 

71 
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Surface Structure of Vitreous Silica - 
A Molecular Dynamics Study 

By S. H. Garofalini 
Department of Ceramics 
Rutgers University 

The importance of vitreous silica ( v-SIOg ) in catalysis, microelectronics, 
optical fibers, and commercial glasses is well appreciated by physical scientists 
working in these areas. Many of the details of the structure and dynamics of 
v-Si0 2 at an atomistic level must be inferred from experimental techniques 
which provide only average properties (e.g. X-ray and neutron diffraction, 

IR and Raman spectroscopy). Also, the amorphous and insu'lative nature of 
v-SiOg have limited the use of the more common surface analytical techniques 
for studying the silica surface in an unambiguous manner. 

The lack of these details at the atomic level provides impetus for the 
application of the molecular dynamics (MD) computer simulation technique to 
study v-Si0 2 . Once the accuracy of this technique has been established, it 
can be used to provide significant information into the structural and dynamic 
behavior of atoms in v-Si0 2 and at the v-Si0 2 surface under a variety of condi- 
tions. 

The excellent computational facilities at NASA-Ames have enabled me to 
perform this evaluation of the adequacy of the MD technique to accurately sim- 
ulate v-Si0 2 * The work last summer showed that the technique accurately simu- 
lates atomic motion in bulk v-Si0 2 (1) . This summer, the work has been aimed 
at establishing if the surface of v-Si0 2 can be accurately simulated. 

The result of the work this summer clearly indicates the adequacy of using 
the MD technique to study the v-Si0 2 surface. The simulations reproduced the 
expected predominance of the oxygen atoms (rather than silicon atoms) at 
the surface, as well as the occurrence of several non-bridging oxygen. An 
appropriate decrease in the S i -0-Si bond angle was also found to occur at the 
surface. The radial distribution function (RDF) as a function of distance from 
the surface was evaluated so as to determine the structural changes occurring at 
the surface. Figure 1 shows the RDF from (a) atoms in the interior of the sample 
{approximating the bulk), (b) all atoms within loR of the surface, and (c) all 


atoms within 5 & of the surface. The first peak in each RDF is due to the 

O 

S1-0 pairs and occurs at the appropriate bond distance of 1.62A. The shoulder 
In the Si-0 peak in Figure lc was found to be due to the several non-bridging 
oxygen (NBO) at the surface, (see Figure 2). 

This decrease in the Si-0 distance between the NBO and Si as compared to 
the Si-0 spacing for the normal bridging oxygen is very Important. Previous 

m 

calculations have Indicated that the Si-0 distance for NBO should be -~0.08A 
less than the normal Si-0 distanced). Our simulation gives this result. 

This result is significant in that it shows that the MD simulation used 
here accurately reproduces surface and defect properties of v-SiOg. 

In addition, a study of the motion of K atoms at the KgO'SSiOg glass sur- 
face immediately after fracture was simulated. Experimental studies (ion 
scattering spectroscopy) have indicated a higher than expected concentration of 
K atoms occurs on the fracture surface. Results of my simulation show that K 
atoms preferentially jump to the fracture surface within a few picoseconds. 

References : 

1. Garofal ini , S. H., J. Chem. Phys. 76(1982) 3189. 

2. Lucovsky, G., Phil. Mag. B, 39(1979) 513. 
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Aerosol Climate Effect Program-Data Comparison 


ORIGINAL PAGE 13 Jindra Goodman 

n c poor QUALITY San Jose State University 

Professor of Meteorology 


The first project was associated with the Aerosol Climate Effects (ACE) 
program. It is well known that stratospheric aerosols play an important 
role in global climate. To assess their effect on the energy budget we have 
to understand their physical and chemical properties, time and space varia- 
tion and evaluate the effects of major volcanic eruptions. Even though 
information on aerosol properties have been gathered over past two decades 
there remains an urgent need for simultaneous measurement using different 
sampling techniques. Therefore as part of the ACE program a continuous 
effort has been made to intercompare instruments measuring the same or 
similar aerosol properties. The sampling platform is the U-2 aircraft and 
the instruments sampling during ACE missions are listed in Table 1 along 
with their characteristics. 

Our specific research interest concentrated on the comparison of the 
Ames wire impactor (AWI) with instruments, which measure similar aerosol 
properties (size distribution). Specifically the AWI was compared with the 
ASAS probe and to a lesser extent with the QCM cascade impacter. The major 
physical difference between the AWI and the ASAS probe is their location. 
The AWI is located directly in the air stream, while ASAS probe must bring a 
sample inside probe. To insure the representation of the sample, several 
conditions have to be fulfiled. Sampling should be isokinetic and isoaxial 
(Fuks, 1975). Detailed analyses of this condition in connection with ASAS 
probe indicates that there may be siime problems. The original mounting 
hindered isokinetic flow, but recently a diffuser cone was installed to 
insure isokinetic conditions. The probe location ori the U-2 doesn't allow 
isoaxial sampling. There is a very sharp bend (^90°) in the intake probe. 
This results in the loss of larger particles and the resulting errors become 
altitude dependent (Stokes number of particle increases with altitude). The 
comparison of the ASAS data with the wire impactor indicated such a trend 
(see Fig. 1). Additional errors in the size determination can be connected 
with variations in the refractive index of particles. On the other hand, 
sampling with the AWI depend on the collection efficiency of the wire (Lem 
and Far low, 1979), sampling volume and the analytical procedure to determine 
particle size. 

Detailed examination of data from the ASAS probe and AWI yield relative- 
ly good agreement, if we take into consideration that a 10% uncertainty in 
particle size will yield a 30-50% error in concentration. 
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Aerosol mass data obtained using the QCM cascade Impactor didn't agree 
well with the AWI. This may be due to the fact that the size intervals are 
larger and the particle size is determined by a 50% collection efficiency. 
Summary of the data are presented in Table 2. 

The second research activity involed the preparation of a proposal to 
study the neutralization of tropospheric acid aerosols by ammonium and their 
participation in the nucleation processes and connection with acid precipita- 
tion. 


References: 


N. A.Fuks, 1975: Review papers-Sampl ing of Aerosols. Atmos. Environment, 

vol . 9, p. 697-707. 

M. Y. Lem and N. H. Farlow, 1979: Efficiency of Aerosol Collection on 
Wires Exposed in the Stratosphere. NASA 81147, 28 pp. 
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TABLE 1 


INSTRUMENT 

ACTIVE SCATTERING 
AEROSOL PROBE-ASAS 

CONDENSATION 
NUCLEI COUNTER 

AMES HIRE IHPACTOR 
AUI 

QUARTZ CRYSTAL 
MICRO BALALNCE (QCM) 

FILTER 

quant ity 
Manured 

aerosol concentration 
In 32 size ranges frc* 
. l^n dia. to 3^m dfa. 

concentration 

concentration 
and size 
distribution 

aerosol mass In 10 
size ranges .OS 1 * 
to 2!ju* 

toL’al S0 4 " aass 


TABLE 2 

1982 ACE FLIGHTS ' 


DATE 

LOCATION 

ALTITUDE 

KM 

No/cm 3 

>0.06vm 

AEPSOL 

VOLUME 

fH 3 /CC 

S0 4 MASS 

m- 

COMPARISION OF SIZE DISTRIBUTION 
OBTAINED FROM ASAS WITH AWI 

COMPARISON 
BETWEEN QCM 
AND AWI 

15 APRIL 

BAJA CA 

16.6 

18.3 

0.1667 

1.0819 

ASAS Missed some large particles 

no data avail, 
yet 

15 APRIL 

BAJA CA 

18.3 

28.9 

.2587 

1.6789 

agree 

. ' 

t 

15 APRIL 

| BAJA CA 

20.7 

6.1 

.2656 

1.7237 

no data 

* 

16 APRIL 

CANADA 

16.8 

18.8 

.1924 

1.2487 

ASAS missed some large particles 

• 

16 APRIL ] 

CANADA 

18.3 

5.1 

.0860 

.5581 

ASAS missed some targe particles 

. i 

16 APRIL 

CANADA 

20.7 

2.2 

.1167 

.7574 

■a data 

m 

19 APRIL 

BAJA CA 

20.7 

2.5 

.071B 

.4659 

no data 

■ 

19 APRIL 

BAJA CA 

20.7 

4.85 

.0954 

.6191 

no data 

• 

5 MAY 

BAJA CA 

15.3 

23.4 

.0603 

.3913 

ASAS reads more small particles 
agreement on large particles 

disagree 

5 MAY 

BAJA CA 

20.7 

3.6 

.0340 

.2206 

agree 

disagree 

7 MAY 

ROC<Y MT. 

16.8 

23.7 

.3666 

2.3792 

ASAS missed some large particles. 

N/A 

7 MAY 

ROCKY MT. 

18.3 

12.6 

.1718 

1.1149 

ASAS missed some large particles 

agree 
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LIFE TESTING OF ACCUMULATOR BLADDERS (AIBS SYSTEM) 

JOHN S. HILTEN 
AUGUST 23, 1982 

EMBRY-RIDDLE AERONAUTICAL UNIVERSITY 
INSTRUCTOR - AERONAUTICAL ENGINEERING DEPARTMENT 


Two Rotor System Research Aircraft (RSRA) have been built by Sikorsky for NASA. 
Each of these aircraft is unique in the sense that there are three flight con- 
figurations: (1) as a helicopter, (2) as an airplane (with the addition of en- 
gines, wings, and a horizontal tail, and with the removal of the rotor blades), 
and (3) as a compound (with the rotor blades back on). In the compound config- 
uration the aircraft flies like a helicopter or like an airplane at the dis- 
cretion of the pilot. 

The basic purpose of the aircraft is to serve as a test vehicle for experimen- 
tal rotor systems. 

Helicopters in general are subject to a wide variety of vibrations; the most 
severe of these is known as the blade passing frequency vibration. The RSRA in 
its present configuration with five rotor blades that typically rotate at about 
200 RPM generates a severe vibration at about 17 HZ. To reduce the transmis- 
sion of this vibration to the rest of the helicopter the transmission is sup- 
ported laterally by 4 active isolators (shock absorbers) that attenuate the 
vibration level. Each isolator has accumulators with bladders that separate a 
pneumatic section from a hydraulic section. 

This research project will take 3 different bladders and life cycle them at 
pressure conditions normally encountered in flight for 12 x 10 6 cycles each; 
this is to verify their durability or mode of failure. 

The envisioned procedure is to pressurize the pneumatic side of the bladder to 
1100 psi with room temperature nitrogen and then to pressurize and cycle the 
hydraulic side via a controller and servo valve to 1650 psi ± 80 psi, at 17 
HZ. The input to the controller will be from a sine wave generator set to 17 
HZ and from the feedback pressure transducer measuring the hydraulic oil 
pressure. A ratio potentiometer will be used to match the two input signals, 
A bias potentiometer will be used to set the static pressure level and a gain 
potentiometer will be used to set the level of the pressure modulation. A 
dither contrdl and limit control will also be provided. The output of the 
controller will be fed into the input of the torque motor which will mechanic- 
ally position the pilot spool valve in the servo valve. A 3000 psi hydraulic 
pump with a maximum flow rate of 31 gal, /min. will furnish the pressure for 
the servo valve. The controller, torque motor, servo valve, and hydraulic 
pump are available commercially. * 


John S. Hil t.en/Embry-Riddle Aeronautical University 
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It is desirable to shut the system down automatically in the event of a bladder 
failure or any other malfunction in the system. Since the system pressure is 
servo controlled, some failure modes would only show a momentary and hard to 
detect pressure change. Other detection methods such as flow, liquid level, 
and current also do not look promising. The following scheme (figure 1) is 
proposed. Once each hour for a period of three minutes (minute 57, 58, and 59) 
the hydraulic pressure to the bladder is shut down. This is done by using a 
clock-driven cam to open a switch that is in series with the torque motorcoil. 
With the hydraulic pressure removed, the captive nitrogen pressure will return 
to 1100 psi if no nitrogen pressure has been lost due to a leak in the accumu- 
lator or a bladder failure. If a failure has occurred, a reed switch in series 
with the hydraulic pump start- stop relay opens shutting the system down. This 
switch is a normally open reed switch rotated to the 1100 psi position on a 
specially modified Bourdon tube pressure gage provided with a magnetic pointer. 
This switch is by-passed by a second clock-driven cam-activated switch for all 
but minute 58 of the hourly cycle. It should be noted that the temperature 
generated by cycling will probably require a pressure set point higher than 
1100 psi . 

A second reed switch rotated to 1650 psi on the Bourdon tube pressure gage is 
also in series with the hydraulic pump start-stop relay. If the pressure 
deviates from the set point of 1650 psi due to power supply, transducer, signal 
generator, controller, torque motor, servo valve, etc. failure the system will 
be shut down. This second reed switch is by-passed by a third clock-driven 
cam-operated switch for minute 56, 57, 58, 59 and 60 of the hourly cycle. 

Thus for the first 55 minutes of each hourly cycle the oil pressure side of the 
bladder is monitored and in the last 5 minutes of each hourly cycle the nitro- 
gen pressure is monitored. If a bladder failure should occur, its time of 
failure may be known no better than 1 hour which represents 1/2 % of the test 
period. The output of the strain gage pressure transducer will also be record- 
ed on a slow strip chart recorder. 

Other details are shown in figure 1. 

At this writing all components have been ordered but delivery schedules are 
typically quoted at 8-12 weeks. 
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Signal Detection Algorithms in the 
Search for Extraterrestrial Intelligence 
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Steven Lord 

Department of Physics and Astronomy 
University of Massachusetts at Amherst 


In efforts to detect radio transmissions from unknown intelligent extra- 
terrestrial sources, any effective search depends on the employment of sensitive 
data processing algorithms. In conjunction with the Ames SETI Group Science 
Team, I have have explored the possibility of using several previously suggested 
techniques (Drake, '8l; Lord, ’8l) and have suggested some new ones as well. 

One established technique, called pseudobinning, is an examination of the 
output of the Multi-Channel Spectrum Analyser at successively coarser resolu- 
tions in hopes of Hatching the bandwidth of any incomming signal. This summer I 
assisted by designing digital circuits for this technique, and by exploring 
theoretical threshhold setting calculations (Drake, ibid.) that prescribe tow to 
achieve a desired false alarm rate using this method under noise conditions. In 
addition, I developed software algorithms for off-line use of this method, and 
determined the assembly language version's minimal time budget. 

Ih the two summers I have provided the group with ideas for various alter- 
native algorithms and processing techniques including the block diagrams for a 
data &y;h&y (time, frequency) transposition scheme that arranges MCSA data in a 
mOff e&cessible form, followed by cross-correlation processor that serves as a 
matched filter for on-again, off -again carrier waves. As a "nonspecific signal" 
type detector, a pattern recognition simulation was written which isolated and 
highlighted detected signals in noise on a video graphics display. The tech- 
nique used an adjancency algorithm, which detected the presence of multiple 
detections closely associated in frequency and/or time. 

To prepare for the day when the available detection equipment includes 
a radiotelescope with dual othogonally polarized receivers, algorithms 
possessing greatly enhanced sensitivity to polarized signals have been developed 
(Lord, ibid.), which, by processing the complex output of two MSSAs, c^namically 
match to the polarization of the incomming signal. The eventual desirability of 
such coherent detection was discussed and agreed upon by the team, and this 
technique called generalized coherence, is likely to be the subject of further 
investigation. 

By attempting to determine the sensitivity of various algorithms, espe- 
cially pseudobinning and coherent detection, I studied the work of Oliver in 
Drake, ibid. This led me to conclude that there did in fact exist a general way 
of determining' unambigously and objectively the sensitivity of an algorithm. 


This method is shovr. in the fig-are belov. Once the algorithm’s response to a 
knovn noise spectrum in learned, either by theory or simulation, as veil as the 
algorithm’s response to relatively large signal to noise ratio environment, a 
single sensitivity estin&te here called S/N may be obtained for any desired pro- 
bability of false alarm and probability of missed signal, here represented by Pr, 

Ti »*• 


This method can be used to rate two candidate algorithms that may have been 
designed to respond to the same signal type. 
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Figure 1. 



Step 1. Use fixed Pr^ to give T. Here 
ve plot 1 minus the cumulative probability 
density faction of the algorithm’s response 
to noise alone. A desired probability of 
a false alarm yields the threshhold T 
of the algorithm (the output level at vhich 
the algorithm vill declare the presence 
of a signal.) 



Step 2. Ve use T and Pr ws to find S/N. 

Here ve plot 1-CPD for various values of 
S/N input to the algorithm. For a fixed 
Pr^ns there is one S/N curve vhich passes 
through the point (T,Pr mJi ). The S/N^. value 
may be thought of as the algorithms sensit- 
ivity for the given Pr^ and Preconditions. 
A lover S/N corresponds to a higher sens- 
itivity. 
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Infrared Spectroscopic Studies of the Effect of Elevated 
* Temperature on the Association of Amino Acids with Clays 


John W. Macklin 

* University of Washington 

Department of Chemistry 

One of the models associated with origin of life chemistry at Ames involves 
the enhancement of peptide bond formation between amino acids by heating in 
the presence of clay minerals. Clay and amino acid mixtures are cycled 
repeatedly through wet cold and dry hot conditions to achieve increases in ^ ^ 
oligomer formation of 40 to 200 times that obtained in the absence of clay. ’ 
The explicit chemical mechanism is not known. Recent chemical studies suggest 
that a covalent type association between the clay and amino acid may constitute 
tiie "activated intermediate" that leads to oligomerization.-^ The research 
described in this report seeks to determine the nature of clay amino acid 
association at elevated temperatures and elucidate the mechanism for pro- 
motion of amino acid oligomerization in the presence of clays. 

The experimental approach is to use Fourier transform infrared spectroscopy 
(FUR) to measure variations in the vibration spectra of clay-amino acid 
samples with increasing temperature. Thin, free standing clay films are 
prepared by evaporation of water from large droplets containing the clay or 
clay-acid mixture placed on a smoothed sheet of flexible polyethylene . 

Silica and alumina samples are measured as fi3.ms on AgCl. Pyroglutamic 
acid was chosen for these studies because the carboxilate group will have 
the same character as that of other amino acids but it is not likely to 
form intermol ocular peptide bonds that would displace covalent associations 
of the carboxilate group with the clay. Hie clay acid association is thus 
preserved for observation. The samples are held undisturbed in a heated 
thermostated cell^ while spectra are measured at temperatures near 24, 60, 

120, 160 and 220 6 C. The temperature is measured with a chrome-alomel 
thermocouple near the sample. Spectra taken at lower temperatures are 
computer substrncted • from higher temperature measurements to obtain the 
explicit changes dependant upon increasing the temperature. 

Results from one of the experiments are shown in Figure 1. The known 
temperature dependent variations in the separate spectra of the clayS and 
pyroglutamic acid 6 arc accounted for so that variations dependent upon 
interactions of pyroglutamic acid with the clay can be determined. The 
experiment has been carried out with several types of clay as well as silica 
gel alurnnia in the form of Gibsite. 

A number of intensity variations are consistently observed in the difference 
spectra of clay-pyroglutamic acid mixtures. For example, note the intensity 


increases at about 1200, 1320 and 1780 cm"- 1 - and decreases near 700, 1245, 

1442 and 1680 cm"'* in the measurements shorn in Figure 1. 'Hie observed 
changes are consistent with the formation of an ester type association 
between the pryoglutamic acid carboxyl group and the mineral surface at 
elevated temperatures. For example, formation of an aceyl silicate may be 
indicated. We are continuing experiments in order to fully define the meaning 
of the spectral variations. 

My experience as a NASA -A, SHF, has been most rewarding. I intend to continue 
ray association with NASA colleagues and am formally proposing to contribute 
to other NASA projects by contractual arrangement. 


‘ References 

1. N. Lahav, D. White and S. Chang, Science , 201 , 67 (1978) 

2. D. White, *.J. Erickson, J. Mol. Evol , , 16 , 279 (1980) 

3. (a) D. Jewett and J. Lawless, Naturwissenchaftcn , 68 , 570 (1981) 

(13) D. White, J. Erickson, private communication, manuscript in 
preparation. 

4. T. Wydcven and M. Leban, Anal . Chcm . 39 , 1673 (1976) 

5. "The Infrared Spectra of Minerals", V. C. Farmer, Chap. 5, pp. 331 -64, 
Mineralogical Society, London (1974) 

6. Marino Vitores ■'Lozano and Anne -Marie Bellocq, J. Chimie Physique, 70 , 

1337 (1973 


ORIGINAL PASS 19\ 
OF POOR QUALITY 


0 2 16 


ORIGINAL PAGE IS 
OF POOR QUALfTY 


Hi 


2-35'-/2'5°C 



2,35 - 2 25“C 



a^5~ /3o*^ 


WRVENUMBEflS CtfrJ. 

TTIR difference spectra of sodium moritmorillonite containing 
50% pyroglutamic acid by weight taken at various temperatures 
The temperatures at which the spectra included in each sub- 
traction were measured are indicated below each trace. 





TITAN'S AEROSOLS 
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Peter H. Smith 

Lunar and Planetary Lab 
University of Arizona 


Most of the observations of Titan that are now available really "see" 
the upper haze layers which are known to blanket the entire lower atmosphere 
and surface. Therefore, when we study Titan with optical methods we are 
learning about aerosol properties. Although it has been three years since the 
first flyby mission to Titan, there is still no general concensus as to what 
the basic properties of the main haze layer, i.e. size, optical constants 
and pressure levels, are. It was in order to reconcile some of the conflicting 
results proposed by different groups working in this area that I carnet to Ames 
and decided to work with James Pollack. 

The observed characteristics of the haze layer are now accepted as follows. 
Ground-based observations of the geometric albedo reveal an increasing amount 
of absorption of solar energy toward shorter wavelengths in the blue. This 
has interpreted as an absorbing haze high in the atmosphere. Voyager I 
confirmed this hypothesis and placed the haze layer about 250 km above the 
solid surface. Both the Pioneer and Voyager missions found this haze to be 
extremely highly polarizing, Tomasko and Smith (1982) interpret this to be 
evidence that the aerosols are small (<''-0.1 /Am), dark and that their size 
increases with depth. Physical models of Titan's aerosols (Toon et ah, 1980) 
also predict this behavior but with larger sizes. Voyager I found a new property: 
the haze is very forward scattering. Rages etal_. (1982) explain that the 
observed forward scattering is indicative of the diffraction from particles of 
radius greater than 0. 1 9y.‘,m and they favor a size nearer 0.5 ytm. Voyager II 
confirms this forward scattering character of the atmosphere leaving investi- 
gators in tho position of choosing between mutually exclusive models for the 
haze layer. 

My first project this summer at Ames has been to put together the Pioneer 
and Voyager data sets in such a way as to derive the single scattering 
properties of an aerosol particle. Normally, these properties serve as the 
input to an atmosphere model. From a preliminary analysis I've found that the 
phase function is decidedly non-spherical : forward scattering yet with little 

back peak and highly polarizing (see Figure). We investigated prolate spheroids 
last summer (Smith, 1982) with some success. Of course, because other shapes 
may be more appropriate, we have started a literature search of the various 
laboratory data pertaining to non-spheres of different classes. Analytical cal- 
culations are only possible for spheroids and cylinders. 

In analyzing the largest phase angle images this summer, it has become clear 
that plane parallel atmosphere calculations are inadequate. So to really define 
the single scattering properties at the most forward scattering geometries, I 
am collaborating with Kathy Rages who has written a spherical scattering program 
and can help develop correction factors to boost the intensities calculated from 
my plane parallel program. We are very close at the end of this summer to a 
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final definition to the forward peak of the phase function for the atmosphere of 
haze aerosols and are ready to begin the search phase of our project. I plan 
to continue this search at the Lunar and Planetary Lab in Tucson and will 
stay in communication with Pollack and Rages as the project progresses. 

In developing the procedures and techniques for analyzing the Titan data 
it has been useful to exchange ideas with other workers in the field here at 
Ames. I'm especially indebted to Brian Toon and Morris Podolak, on leave from 
Tel Aviv. 
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